Volume 104, Number 1, January-February 1999 

Journal of Research of the National Institute of Standards and Technology 

[J. Res. Natl. Inst. Stand. Technol. 104, 11 (1999)] 

Thermodynamic Temperatures of the Triple 

Points of Mercury and Gallium and in the 

Interval 217 K to 303 K 



Volume 104 



Number 1 



January-Febraary 1999 



M. R. Moldover, S. J. Boyes', C. W. 
Meyer, and A. R. H. Goodwin^ 

National Institute of Standards and 

Technology, 

Gaithersburg, MD 20899-0001 



We measured the acoustic resonance fre- 
quencies of an argon-filled spherical cav- 
ity and the microwave resonance frequen- 
cies of the same cavity when evacuated. 
The microwave data were used to deduce 
the thermal expansion of the cavity and 
the acoustic data were fitted to a tempera- 
ture-pressure surface to deduce zero- 
pressure speed-of-sound ratios. The ratios 
determine {T—T^), the difference be- 
tween the Kelvin thermodynamic tempera- 
ture T and the temperature on the Inter- 
national Temperature Scale of 1990 
(ITS-90). The acoustic data fall on six 
isotherms: 217.0950 K, 234.3156 K, 
253.1500 K, 273.1600 K, 293.1300 K, 
and 302.9166 K and the standard uncertain- 
ties of (T—Tijo) average 0.6 mK, depend- 
ing mostly upon the model fitted to the 
acoustic data. Without reference to 
ITS-90, the data redetermine the triple 
point of gallium Tg and the mercury 
point r„ with the results: T^/T^ = 
(1.108 951 6 ± 0.000 002 6) and TJT„ = 
(0.857 785 5 ± 0.000 002 0), where 
r„ = 273.16 K exactly. (All uncertainties 
are expressed as standard uncertainties.) 
The resonator was the same one that had 



been used to redetermine both the universal 
gas constant R , and T„. However, the 
present value of Tg is (4.3 ± 0.8) mK larger 
than that reported earlier. We suggest 
that the earlier redetermination of Tg was 
erroneous because a virtual leak within 
the resonator contaminated the argon used 
at Tg in that work. This suggestion is 
supported by new acoustic data taken when 
the resonator was filled with xenon. For- 
tunately, the virtual leak did not affect the 
redetermination of R . The present work 
results in many suggestions for improving 
primary acoustic thermometry to achieve 
sub-millikelvin uncertainties over a wide 
temperature range. 
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1. Introduction 



In the Introduction, we briefly review the historical 
context of the present work, the conceptual basis of 
primary acoustic gas thermometry with spherical res- 
onators, the present results, and the significant compo- 
nents of the uncertainty of these results. 
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1.1 Historical Context 

Here we report new values for the thermodynamic 
temperatures of the triple points of mercury Tm and 
gallium Tg and the difference (T—T^) between the 
Kelvin thermodynamic temperature scale T and the In- 
ternational Temperature Scale of 1990 (ITS-90) from 
217 K to 303 K. This work documents significant pro- 
gress in a program at NIST/NBS to exploit spherical 
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cavities for primary acoustic gas thermometry. This pro- 
gram began during 1978 when Moldover et al. [1] mea- 
sured the frequencies of the acoustic resonances in an 
argon-filled spherical cavity and also deduced the ra- 
dius of the cavity from the frequencies of microwave 
resonances within it. In doing so, they demonstrated the 
essential elements of primary acoustic thermometry us- 
ing a spherical cavity. Important advances were made 
by Mehl and Moldover [2] and by Moldover, Mehl, and 
Greenspan [3] who published a detailed theory for the 
acoustic resonances of a nearly-spherical, gas-filled cav- 
ity as well as extensive experimental tests of the theory. 
These results guided Moldover et al. [4] in assembling a 
3L, steel-walled, spherical cavity sealed with wax (the 
"gas-constant resonator") which they used during 1986 
to redetermine the universal gas constant R with a rela- 
tive standard uncertainty of 1.7 X 10"'', a factor of 5 
smaller than the uncertainty of the best previous mea- 
surement. Mehl and Moldover [5] also developed the 
theory of nearly-degenerate microwave resonances in a 
nearly-spherical cavity and showed how to use a few 
microwave resonances to deduce the volume of the cav- 
ity. Their theory was tested by Ewing et. al [6] who 
showed that a microwave measurement of the thermal 
expansion of the gas-constant resonator from 273 K to 
303 K was consistent with a measurement based on 
mercury dilatometry. 



The gas-constant resonator had not been optimized 
for the determination of the thermodynamic tempera- 
ture T; however, it was very well characterized and it 
was available for measurements prior to the scheduled 
replacement of the International Practical Temperature 
Scale of 1968 (IPTS-68) with ITS-90. Thus, Moldover 
and Trusler [7] used the gas-constant resonator during 
1986 to determine Tg, the thermodynamic temperature 
of the triple point of gallium. Subsequently, during 
1989, two of the present authors (M.R.M. and C.W.M.) 
used the gas-constant resonator to study the temperature 
scale in the range 213 K to 303 K and found evidence 
that the redetermination of Tg reported by Moldover and 
Trusler [7] was in error by approximately 4 mK [8]. 

Here, we report the results of more recent measure- 
ments conducted with the gas-constant resonator, pri- 
marily during 1992, which lead to new values of 
(T-Tgo) on the five isotherms: 217.0950 K, 234.3 156 K, 
253.1500 K, 293.1300 K, 302.9166 K. One of these 
isotherms (302.9166 K) is very near Tg and another 
(234.3156 K) is very near T„„ the thermodynamic tem- 
perature of the triple point of mercury. For the present 
data, the standard uncertainty of (T—Tgo) is approxi- 
mately 0.6 mK, depending mostly upon the model fitted 
to the acoustic data. This sub-millikelvin uncertainty is 
much smaller than the uncertainties characteristic of 
other methods of primary thermometry in this tempera- 
ture range. (See Fig. 1 .) 
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Fig. 1. Comparison of the Kelvin thermodynamic temperature scale and ITS-90. The present 
data are designated 'This Work (1994)" and the solid curve is a spline fit to them. The other 
data sources are: NBS gas thermometry, Ref. [15]; NML gas thermometry Ref. [12]; corrected 
PRMI gas thermometry, Ref. [13]; NPL total radiation thermometry, Ref. [14]; NBS acoustic 
thermometry, Ref. [7]; UCL acoustic thermometry, Ref. [11]. 
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The present value of T^ is (4.3 ± 0.8) mK larger than 
the Moldover-Trusler [7] value. (All uncertainties in this 
manuscript are expressed as standard uncertainties.) 
With the benefit of both additional data and hindsight, 
we conjecture that the Moldover-Trusler determination 
of Tg was erroneous because the argon used in their 
work was progressively contaminated, perhaps by vir- 
tual leak into the resonator. (Any volume that was sealed 
from the laboratory and was connected to the resonator 
by a path of low pumping speed would act as a virtual 
leak. If such a volume were exposed to a contaminating 
gas (e.g., air at 100 kPa), it would fill rapidly via 
Poiseuille flow. Subsequently, the volume would be a 
persistent source of contamination because it would take 
a very long time to empty via molecular flow at low 
pressure.) Below, we present some data acquired when 
the gas-constant resonator was filled with xenon that 
strongly support this conjecture. The need to conjecture 
calls attention to a significant weakness of the gas-con- 
stant resonator and the apparatus associated with it: 
there were no satisfactory provisions for detecting con- 
tamination of the thermometric gas after it had been 
admitted into the resonator. Fortunately, all of the re- 
sults from the gas-constant resonator on the 273.16 K 
isotherm are mutually consistent; thus, there is no evi- 
dence that contamination was a problem during the re- 
determination of R . 

At the conclusion of the present work, the gas-con- 
stant resonator was disassembled. The dimensions of its 
component hemispheres were re-measured. Now the 
hemispheres and the associated apparatus are being re- 
constructed to correct deficiencies uncovered in the 
present work and to extend the range of primary acous- 
tic gas thermometry up to 750 K. Improvements to the 
apparatus and procedures include: (1) flowing the ther- 
mometric gas through the apparatus to reduce the resi- 
dence time of the gas by a factor of 1 00, thereby reduc- 
ing exposure to possible contamination; (2) tuning the 
ports that admit gas to the resonator [9] thereby attain- 
ing simultaneously a high acoustic impedance near ra- 
dial resonances and a low flow impedance; (3) bakeable 
gas-handling system and transducers (no polymer seals) 
minimizing possible contamination; (4) analysis of the 
gas exiting the resonator via gas chromatography; (5) 
simultaneous measurements of acoustic and microwave 
resonance frequencies; (6) positioning the microwave 
coupling probes to optimize the resolution of the nearly- 
degenerate microwave resonance frequencies; (7) provi- 
sions for measuring possible horizontal temperature 
gradients in the resonator; (8) provisions for measure- 
ment of the resonator's temperature on ITS-90 with up 
to five long-stemmed standard platinum resistance ther- 
mometers that can be conveniently inserted and re- 
moved from contact with the resonator; and (9) extended 



data runs to acquire data on many, relatively closely- 
spaced isotherms. (The gas constant apparatus had to be 
disassembled to calibrate the thermometers. Thus, ex- 
tended data runs were precluded by the concern that the 
thermometers might drift between calibrations.) This 
list of projected improvements makes it clear that the 
present results will not comprise NIST's final contribu- 
tion to acoustic gas thermometry in the temperature 
range 217 K to 303 K or above. 

1.2 Conceptual Basis of the Present Measurements 

1.2.1 Acoustic Gas Thermometry 

Primary acoustic thermometry relies on the connec- 
tion between the speed of sound in a gas and the thermo- 
dynamic temperature. Elementary considerations of hy- 
drodynamics and the kinetic theory of dilute gases leads 
to relationships between the thermodynamic tempera- 
ture T, the average kinetic energy E in one degree of 
freedom, and the speed of sound u : 



3£ = 



1 



mVrms ■ 



:kT, 






(1) 



Here, Vrms is the root mean square speed of a gas 
molecule, m is its mass, k is the Boltzmann constant, and 
y is the ratio of the constant pressure to constant volume 
specific heat capacities which is exactly 5/3 for perfect 
monatomic gases. The International System of Units 
assigns the exact value 273.16 K to the temperature of 
the triple point of water T„. From this assignment and 
from Eqs. (1), the Kelvin thermodynamic temperature T 
of a gas can be determined from the zero-pressure limit 
of the ratio of speed of sound measurements at T and T„ 
with the equation 



273.16 K 



lim ( u^(p,T) 



(2) 



In principle, Eq. (2) could be used to calibrate ther- 
mometers on the Kelvin thermodynamic scale and ulti- 
mately, this may become accepted practice. 

Today, some useful implementations of Eq. (2) can be 
accomplished with stable thermometers that have never 
been calibrated using ITS-90. For example, in this work 
we could have used uncalibrated ("transfer-standard") 
thermometers to redetermine the thermodynamic tem- 
perature of the triple point of mercury T^ (and of other 
fixed points). To do so, we would have recorded the 
indications of the uncalibrated thermometer when it was 
in a mercury-point cell and when it was at r„. Subse- 
quently, the thermometer's indications would be used to 
help adjust the temperature of the thermometric gas to 



13 



Volume 104, Number 1, January-February 1999 

Journal of Research of the National Institute of Standards and Technology 



Tn, and then to T„ where u^ip ,T„,) and u^{p ,T„) would be 
measured. Finally Eq. (2) would be applied. In practice 
we used standard platinum resistance thermometers 
(SPRTs) that had been calibrated on ITS-90 to adjust 
the temperature of the argon. We used SPRTs for two 
reasons. First, SPRTs are known to be very stable and 
second, our results at temperatures between fixed points 
could be preserved and reported as the difference 
{T— Tgo) as a function of T or of Tgo- 

1.2.2 Acoustic and Microwave Resonances in a 
Spherical Cavity 

In this work, the speeds of sound in argon and in 
xenon were deduced from measurements of F^, the fre- 
quencies of the radially-symmetric acoustic modes in 
the gas-constant resonator. (Conventionally, F^ =f^ + ig^ 
is a complex number such that /a is the center frequency 
of the resonance and g^ is the half-width of the reso- 
nance. In the context of frequency measurements, the 
subscripts "a" and "m" are used to distinguish acoustic 
modes from microwave modes; in the context of triple- 
point temperatures, the subscripts "a" and "m" identify 
argon and mercury, respectively.) There is a well-devel- 
oped theory for the radially-symmetric acoustic modes 
that has been confirmed by detailed experiments [2], 
[4], [10]. The frequencies of these modes are only 
weakly sensitive to deformations of the shape of the 
cavity from a perfect sphere so long as the volume of 
cavity V{p,T) is unchanged. Thus, accurate measure- 
ments of u^{p,T) do not require accurate measurements 
of the shape of the cavity. They do require an accurate 
measurement of V{p,T), which is usually much easier. 
(In Ref. [4], the very small pressure dependence of 
V(p ,T) was calculated and checked by a measurement. 
For the remainder of this Introduction, the pressure de- 
pendence of V will be ignored.) We write u{p,T) in 
terms of the frequencies and the volume, 

u{p,T)= \f.(p,T) + ^f.(p,T)] X [ViT)f"/A^. (3) 

In Eq. (3), the term AfJj>,T) represents several small, 
theoretically based, mode dependent corrections which 
must be added to the measured resonance frequencies 
and which depend mostly upon the boundary conditions. 
The constant A^ is known exactly; it is an acoustic eigen- 
value that depends upon the mode for which f^{p ,T) is 
measured. 

As shown in Ref. [5], the microwave modes within a 
nearly-spherical cavity occur in nearly-degenerate mul- 
tiplets with 2/ -n 1 components. (/ is a positive integer.) 
The frequency of each component of a multiplet de- 
pends upon the details of the shape of the cavity; how- 
ever, the average frequency of each multiplet is not 



sensitive to smooth deformations of the cavity that leave 
its volume unchanged. In analogy with Eq. (3), the 
speed of light in the gas c{p,T) can be related to the 
frequencies of the microwave resonances F^ =fm + igm 
through 

c(p,T) = (f^ipj) + AUipJ)) X [V{T)f'VA^. (4) 

In Eq. (4), the term Af^{p,T) represents the mode-de- 
pendent corrections that must be added to the measured 
resonance frequencies, A„, is an eigenvalue for the mi- 
crowave multiplet chosen, and the brackets "(. . .)" de- 
note the average over the components in a multiplet. 

One can combine Eq. (3) and Eq. (4) to eliminate 
V{T) and obtain 



u(p,T) _ f,(p,T) + Af,(p,T) A^ 

c(p,T) (f^(p,T) + Af^(p,T)) A,- ^^^ 



For primary thermometry in an ideal situation, Eq. (5) 
is used at the temperature T and at T^ and the zero-pres- 
sure limit is taken at each temperature. This leads to an 
idealized working equation: 



T lim ( f,(p,T) + Af,(p,T) 

273.16 K p^O\{f^(p,T) + Af^(p,T)) 



X 



lim / {/'n.(p,rj + A/„(p,r„)) 



(6) 



As implied by the absence of V{T) in Eqs. (5) and (6), 
and as emphasized by Mehl and Moldover [5], the spher- 
ical cavity plays a limited role in measuring ulc, the 
ratio of the speed of sound to the speed of light. One may 
view the cavity as a temporary artifact that must remain 
dimensionally stable just long enough to measure /^(p) 
and (fmip)) at the temperature T and that must not 
change its shape (and eigenvalues) too much when the 
frequency measurements are repeated at T^. (Small, 
smooth changes in the shape of the cavity affect the 
eigenvalues only in the second order of the small 
change.) Although Eq. (6) fully accounts for the change 
of the cavity's volume upon going from T to T^, the 
computation of the correction terms AfJj},T) and 
Af^{p,T) require moderately accurate knowledge of the 
cavity's dimensions, the electrical conductivity and the 
mechanical compliance of the shell bounding the cavity, 
and other properties of the both resonator and the ther- 
mometric gas. 

1.2.3 Particulars of This Work 

Unfortunately, the gas-constant resonator did not have 
provision for simultaneous acoustic and microwave mea- 
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surements that would be required to use Eq. (6). There- 
fore, as in Ref. [6], V{T) was deduced from measure- 
ments of/ni(r) made when the cavity was evacuated and 
V{T) was assumed to be a reproducible function of T for 
intervals of weeks. This assumption is supported below 
by the important observation that the values of f^(p ,r„) 
that were obtained during the measurement of R in 
1986, the measurement of Tg in 1987, and subsequent 
thermometry in 1989 and in 1992 (this work) agree 
within 1 part in 10^. 

In principle, u (and T) could be determined from 
measurements of the frequencies of a single acoustic 
mode and a single microwave triplet. In this work, as in 
the redetermination of R , five non-degenerate acoustic 
modes spanning a frequency ratio of 3.8 : 1 were used. 
Also, three microwave triplets spanning a frequency 
ratio of 3.4 : 1 were used. The redundant acoustic and 
microwave measurements were used to determine some 
components of the uncertainty in measuring (T—Tgo) 
and to search for limitations of the theories for the cor- 
rections A/a and A/n,. 

The microwave measurements avoid an assumption 
that is often made in gas thermometry; specifically, that 
the volumetric expansion of an assembled cavity is iden- 
tical with that computed from the measured linear ther- 
mal expansion of other artifacts made from the same 
material. In fact, the present microwave measurements 
provide evidence that the expansion of the gas-constant 
resonator was anisotropic, an effect that would not have 
been detected by measurements of the thermal expan- 
sion of one dimension of an artifact made of the same 
metal. 

The extrapolation of / to zero pressure was accom- 
plished by fitting the acoustic data on all the isotherms 
simultaneously to polynomial functions of the pressure 
with temperature-dependent coefficients. This proce- 
dure together with the practice of measuring/n, when the 
resonator was evacuated and the assumption that the 
dimensions of the gas-constant resonator were stable in 
time leads to a modification of Eq. (6) that more nearly 
represents the working equation used in this work: 



Mm (Mp,T) + AMp,T) 



273.16 K p-^0\f,(p,T^) + Af,(p,T^). 



{f^(T) + Aun j 



(7) 



Typically, a single measurement of/ with the gas-con- 
stant resonator had a repeatability of 0.2 X 10"^ /a cor- 
responding to 0. 1 mK at 300 K. The repeatability mea- 
surement of /„ was less than 0. 1 X 10~^/rn. The excellent 
repeatability of these frequency measurements indicates 
that acoustic gas thermometry has the potential to be 



very accurate. The larger uncertainties encountered in 
this work are listed in Sec. 1.4 and their evaluation is 
discussed throughout the body of this manuscript. 

1.3 Results 

1.3.1 The Differences (T-Tg^) 

The present results for (T—Tgo) are listed in Table 1 
and they are compared with other results in the same 
temperature range in Fig. 1. Included in Fig. 1 are re- 
sults from primary acoustic thermometry at University 
College London (labeled "UCL Acoustic") at 240 K and 
300 K. The UCL results agree with the present results 
within the remarkably small combined uncertainties. At 
the time of this writing, the UCL results have not been 
published [11]. They are a product of research that 
started after the NIST/NBS program had begun and the 
techniques used at UCL are described in connection 
with a calculation by Ewing et al. [10] of the effects on 
acoustic resonance frequencies of imperfect thermal ac- 
commodation at the shell-gas boundary. 



Table 1. The difference T— Tc 



90 





Isotherm fits 


Surface fit 


Recommended 




linear AAi(7') 






T90/K 


(T-T9o)/mK 


(T-T9o)/mK 


(T-T9o)/mK 


Argon 


302.9166 


3.95±0.73 


4.61 ±0.33 


4.6±0.6 


293.1300 


2.60±0.75 


3.23±0.31 


3.2±0.6 


253.1500 


-2.84±0.60 


-2.29±0.27 


-2.3±0.6 


234.3156 


-3.73±0.60 


-2.92±0.19 


-2.9±0.6 


217.0950 


-4.03±0.67 


-3.55±0.29 


-3.6±0.6 


Xenon" 


302.9166 




4.38±0.66 





' For xenon, the imcertainty is the quadrature sum of 0.40 mK from 
fitting the acoustic data, 0.31 mK from the non-acoustic items in 
Table 2, and 0.43 mK from the virtual leak correction. 

Fig. 1 includes data that were used to establish 
ITS-90. Among these, the present results below the 
triple point of water (T^) agree best with the NML gas 
thermometry [12], and with the corrected PRMI gas 
thermometry [13]. (Before the 1995 correction, the 
PRMI gas thermometry results were much closer to the 
Tgo baseline.) The data on Fig. 1 labeled "NPL Total 
Radiation" came from Ref. [14]; the data labeled "NBS 
Gas" came from Ref. [15]. 

In Fig. 1, the Moldover-Trusler [7] redetermination 
of Tg is labeled "NBS Acoustic, 1988." The difference 
between this value and the present result is 
(4.3 ± 0.8) mK. As mentioned above, we conjecture 
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that the Moldover-Trusler result was erroneous because 
a virtual leak contaminated the argon used in their work. 
In Sec. 8.4 below, we discuss the evidence supporting 
this conjecture and we present some data acquired when 
the gas-constant resonator was filled with xenon that 
also support this conjecture. 

Table 1 summarizes the present results for (T—Tgo) 
that were obtained using argon and xenon. The columns 
labeled "isotherm fits" and "surface fit" resulted from 
two different analyses of the acoustic data for argon. For 
the isotherm fits, the squares of speeds of sound m^ on 
each of the five isotherms listed as well as the isotherm 
273.16 K were fitted by polynomial functions of the 
pressure p. In this analysis, 24 parameters were used to 
represent the six isotherms comprising the u^(p,T) sur- 
face. For the surface fit, three constraints were imposed 
on the u^{p,T) surface, leaving only 11 or 12 parameters 
to be fit by the same data. The constraints applying to 
the temperature range 217 K < T < 303 K were: (1) the 
difference AAi(r) = Ai(r,expt.) — Ai(r,calc.) must be 
an empirical linear function [here, Ai(r,expt.) is the 
coefficient of p in the polynomial and Ai(r,calc.) is the 
prediction from a semi-empirical interatomic potential 
taken from the literature], (2) the coefficient A2(r) of p^ 
in the polynomial is represented as an empirical 
quadratic function of 1/7", and (3) the thermal accom- 
modation coefficient h is exactly one. Constraints (1) 
and (2) are plausible because the present isotherms are 
well above the critical temperature of argon (1.4 < TIT^ 
< 2.0) and well below the critical density (p/pc — 0.02) 
where the virial coefficients of argon are only weakly 
temperature dependent. The third constraint is the tem- 
perature-independent upper bound to h, which will be 
discussed in Sec. 8.3. 

As expected, the constrained (surface) fits resulted in 
smaller statistical uncertainties of the parameters; how- 
ever. Table 1 shows that the constrained and uncon- 
strained fits yield consistent results for (T—Tgo)- Fur- 
thermore, the standard deviations of the fits it{u^) 
relative to (mq^) were comparable. (For the constrained 
fit, 10" o-(m^)/(mo^)= 1.12; for the isotherm fits, 
0.83 < 10" a{uyW) < 1.50.) Thus, the experimental 
evidence does not contradict the constraints. Remark- 
ably, the constrained fit leads to values of (T—Tgo) that 
average 0.6 mK larger than the unconstrained fit. Im- 
posing the constraints increases the zero-pressure limit 
of u^{p,T) on all of the isotherms except the 273.16 K 
isotherm. We cannot explain this; however, in our opin- 
ion, it is advantageous that the constrained fit does not 
give the speed-of-sound data on the 273.16 K isotherm 
a privileged role such that they affect all of the values of 

(^-^90). 



We considered a variety of additional fits (Sec. 8.4) 
including a constrained fit in which AAi{T) was repre- 
sented as a quadratic function of T, and both constrained 
and unconstrained fits excluding acoustic data for the 
(0,6) mode, the mode that deviated most from the aver- 
age of the other modes. These analyses led to values of 
(T—Tgo) that were either between the constrained and 
unconstrained results or very close to them. Obviously, 
the acoustic data could have been analyzed with many 
other combinations of constraints. The data are tabu- 
lated in appendices, allowing the reader to impose the 
constraints that he/she prefers. 

In view of these factors, we recommend the values of 
(T—Tgo) resulting from the constrained fit; however, the 
values are model-dependent. With some arbitrariness, 
we recommend enlarging the standard uncertainty of 
(T—Tgo) to 0.6 mK. The value 0.6 mK encompasses the 
extreme values of (T—Tgo) resulting from the models 
that we investigated and is approximately twice the un- 
certainties resulting from the surface fits. 

We conclude this discussion of the results by noting 
that the recommended uncertainty of (T—Tgo), namely, 
0.6 mK, is very small in the context of results obtained 
by other methods, as shown in Fig. 1 . 

1.3.2 Triple Points of Gallium and Mercury 

The isotherms near the triple points of gallium 7g, 
mercury r^, and water T^ have a special status. For these 
isotherms, our platinum thermometers were used to in- 
dicate when the temperature of the resonator had been 
adjusted to be equal to the temperature of one of these 
triple points. For this very limited function, the ther- 
mometers had to be stable; however, they did not ever 
have to be calibrated on ITS-90. Thus, the present 
acoustic and microwave data redetermine T^ and T^, 
without reference to ITS-90. Using Eq. (7) with the 
microwave data and the constrained fits to the acoustic 
data leads to the results 7^/7;, = (1.108 951 6 ± 
0.000 002 6) and TJT„ = (0.857 785 5 ± 0.000 002 0) 
which determine both 7g and T„, through the definition: 
rw^ 273.16 K, exactly 

1.4 Components of the Uncertainties 

Table 2 lists the important components of the standard 
uncertainty (u^) in the determination of (7"— 7'9o)/rfrom 
the measurements of the quantities in Eq. (7). The eval- 
uation of these contributions is a major portion of the 
body of this manuscript. Here, we outline the phenom- 
ena that contributed to u.. 
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Table 2. Standard uncertainties u,XW^ from various sources in the re-determination of {T—Tgo)IT. The square root of the sum of the squares (RSS) 
is calculated twice: first, including Rows 4 and 5, but not Rows 6 and 7; and second, including Rows 6 and 7, but not Rows 4 and 5 


Source 


217 K 234 K 


253 K 




293 K 






303 K 


Microwave values for [a{T)la{T„)f 


1. Discrepancies among triplets 

2. 5JT) calc. from resistivity (0.04Xp) 

3. 5„(7')calc.-S„(7')meas. 


0.32 0.27 
0.39 0.25 
0.24 0.24 


0.25 
0.12 
0.24 




0.27 
0.12 
0.24 






0.33 
0.18 
0.24 


Acoustic (isotherm fits) 


4. Uncertainty of AoCrwVa^ 

5. Uncertainty of AoCryu^ 


1.67 1.67 
2.31 1.85 


1.67 
1.35 




1.67 
1.67 






1.67 
1.43 



6. Uncertainty of Ao(7'„)/a^ 

7. Uncertainty of AoCryu^ 



Acoustic (surface fit) 

0.25 0.25 

0.52 0.38 



0.25 
0.36 



0.25 
0.31 



0.25 
0.38 



Thermometry 



8. SPRT & bridge repeat. @ T^ (10 (juXi) 

9. Difference between calibrations 

10. Temperature gradient 

1 1 . Non-uniqueness of ITS-90 



0.46 0.43 


0.40 


0.34 


0.33 


0.22 0.15 


0.10 


0.54 


0.78 


0.46 0.43 


0.40 


0.34 


0.33 


0.9 0.0 


0.8 


0.6 


0.0 


Additional sources 


0.20 0.20 


0.20 


0.20 


0.20 


0.21 0.12 


0.05 


0.03 


0.07 



12. Thermal conductivity (0.3 %) 

13. Uncertainty of pressure zero 



14. Isotherm fits: RSS 


3.13 


2.62 


2.40 


2.58 


2.43 


15. Surface fit: RSS 


1.42 


0.92 


1.16 


1.11 


1.13 


16. Isotherm fits: RSSX(7'/mK) 


0.68 


0.61 


0.61 


0.76 


0.74 


17. Surface fit: RSSX(7'/mK) 


0.31 


0.22 


0.29 


0.33 


0.34 



1.4.1 Microwave Measurements 

To determine T, we required the combination of mi- 
crowave frequencies [(fAT^)+^fJTJ)l(f^{T)+^fJT))f 
which equals the ratio [a{T)la{T^y\^ when a{T) is de- 
fined to be the average radius of the spherical cavity 
such that Vo{T) = (4/3)T7[a(r)]' and where Vo{T) is the 
zero-pressure limit of the volume. The ratio [a{T)l 
a{T^y\^ was computed from a polynomial function of 
the temperature that had been fitted to the data for 
(fnff) + A/m(r)). We denote the relative standard uncer- 
tainty of a{T) by Mr(a). The primary components of 
Mr(a) arise from the different values of (fm(T) + A/m(r)) 
for the three microwave triplets studied and from the 
uncertainties in the correction term A/m(r). The differ- 
ences among the triplets account for Row 1 of Table 2. 
Theoretically, the correction term A/n,(r) is proportional 
to 5n,(r), the microwave "penetration depth" or "skin 
depth." Two methods were used to obtain 5^(7"). For the 



first method, 5^(7") was computed from published val- 
ues for the electrical resistivity of stainless steel and 
Row 2 of Table 2 accounts for the uncertainties in adapt- 
ing the dc resistivity data to the present circumstances 
(Sec. 3.2.2.). This approach sets a lower bound to 5n,(r). 
A plausible upper bound for 5n,(r) was calculated from 
the measured half-widths g^ of the microwave multi- 
plets. (The lower-bound values of SJJ') were used in 
computing l^fm{T) for determining the temperature in 
Table 1 ; if the upper-bound had been used, TIT^ would 
have been changed by less than 0.3 X 10"''.) The differ- 
ences between the upper and lower bounds of 8J(T) 
were used to calculate a contribution to the standard 
uncertainty Us{T) which, expressed as a fraction of 
(T—Tgo)/T, appears in Row 3 of Table 2. Negligible 
(< 0.1 X 10 '') contributions to u^a) came from the un- 
certainty of each measurement of /„ and from the uncer- 
tainty of the temperature during the measurement of 
UT). 
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1.4.2 Acoustic Measurements 

The relative uncertainty in the measurement of an 
acoustic resonance frequency was usually less than 
0.3 X 10"''. With limited and well understood excep- 
tions, the five acoustic modes yield results that are con- 
sistent at this level. However, the zero pressure limit of 
the speed of sound was determined by fitting a function 
of pressure to the measured resonance frequencies. The 
correlations among these parameters contribute to the 
uncertainty in the speed of sound ratios accounting for 
the different uncertainties resulting from the isotherm 
fits (Table 2, Rows 4 and 5) and the surface fit (Table 2, 
Rows 6 and 7). 

1.4.3 Temperature Measurements 

In the present work, three capsule-style standard plat- 
inum resistance thermometers were calibrated at the 
triple points of argon, mercury, water, and gallium (T^, 
T^, Tw, Tg) and then installed in the resonator. The un- 
certainty of the thermometry resulted from the uncer- 
tainty of each calibration measurement (Table 2, Row 8) 
and from drifts in the system (thermometers + resistance 
bridge + standard resistor) during the weeks between 
calibrations. The latter was estimated (Table 2, Row 9) 
from the change of the calibrations during the interval in 
which the acoustic measurements were made. Addi- 
tional uncertainty resulted from the small temperature 
difference between the thermometers embedded in the 
top ("north pole") and bottom ("south pole") of the 
spherical shell. The difference never exceeded 0.5 mK. 
We estimate that our imperfect knowledge of the volume 
average of the temperature distribution within the ther- 
mometric gas is no more than 0.1 mK relative to the 
thermometer calibrations. This effect accounts for Row 
10 of Table 2. Except at the calibration temperatures, the 
non-uniqueness of ITS-90 contributed to the uncer- 
tainty of the temperature measurements. (Table 2, Row 

11) 

The "additional sources" of uncertainty in Table 2 
arise from the uncertainty of the thermal conductivity of 
the argon and from imperfect pressure measurements. 
They are considered elsewhere in this manuscript. 

1.5 Organization of This Manuscript 

The remainder of this manuscript is organized into 
major sections as follows: Sec. 2, modifications of the 
apparatus and the experimental procedures; Sec. 3, mi- 
crowave measurements and their interpretation; Sec. 4, 
acoustic measurements; Sec. 5, thermometry; Sec. 6, 
characterization of the gases; Sec. 7, pressure and other 
thermophysical quantities; Sec. 8, analysis of the acous- 



tic results; Sec. 9, comparisons of acoustic results with 
previous results; Sec. 10, other tests for systematic er- 
rors; Sec. 11, tabulated data; and Sec. 12, references. 



2. Apparatus and Procedures 

The apparatus and procedures used have been de- 
scribed in extensive detail in earlier publications. [4,7] 
These details will not be repeated here; instead, we 
discuss the few modifications that were required to im- 
prove the thermostatting for the present work and the 
results of the dimensional measurements that were made 
when the apparatus was disassembled at the conclusion 
of the present work. Figure 2 shows a cross-section of 
the gas-constant resonator and the pressure vessel that 
enclosed it. The pressure vessel was immersed in an 
insulated, thermostatted, well-stirred, methanol bath. 



Clean gas 
and vacuum 




Fig. 2. Schematic cross-section of the gas-constant resonator and 
pressure vessel. 
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2.1 Modifications to the Apparatus 

Several modifications were made to improve the ther- 
mostatting of the gas-constant resonator. When the gas 
constant was measured, a 2 cm long bellows led from the 
isolation valve atop the resonator to the gas handling 
system. (See Figs. 3 and 5 in Ref. [4].) Subsequently, 
during the manipulations associated with the microwave 
measurements, the bellows was damaged. Prior to the 
acoustic measurements of 1989 the bellows was re- 
placed with a 23 cm long, copper tube that had been 
bent into a circle and placed in a horizontal plane nearly 
concentric with the valve atop the resonator. 

During preliminary measurements performed in 
1989, we found that as the temperature of the methanol 
bath was reduced below ambient temperature, the verti- 
cal temperature gradient across the resonator increased. 
When the methanol reached 213 K, the top of the res- 
onator was 3 mK warmer than the bottom suggesting 
that a thermal link from the top to ambient temperature 
existed. This gradient was not reduced by three changes: 
(1) thinning the supports of the pressure vessel, (2) 
improving the radiation shields in the tubes leading to 
the resonator and, (3) improving the stirring of the bath. 

However, the gradient was reduced to about 1 mK by 
surrounding the resonator with a cylindrical heat shield 
comprised of 3 mm thick copper strips. The strips were 
separated from each other but all were thermally an- 
chored to the top and bottom of the resonator with thick 
aluminum strips. The shield was insulated from the 
walls of the pressure vessel by a 3 mm wide space filled 
with argon. This arrangement was used for the mi- 
crowave measurements of the thermal expansion; how- 
ever, the 1 mK gradient at 213 K was larger than desired 
for the acoustic measurements. 

Immediately prior to the measurements of 1 992 three 
additional changes were made to reduce the gradient: 
(1) The copper tube was replaced with a 23 cm long, 
thin-walled, stainless steel bellows (see Fig. 2) which 
was mounted in the same manner as had been the copper 
tube; (2) closed cell foam insulation was glued to the 
underside of the aluminum plate atop of the fluid bath 
to reduce the temperature gradient within the bath itself; 
and (3) a new shield, constructed entirely of copper, was 
installed. This shield, labeled "Isothermal shield" in Fig. 
2, was comprised of a 2 mm thick cylinder and two 5 
mm thick end plates. The shield was suspended from the 
central support tube, but otherwise was not in contact 
with either the resonator or the pressure vessel. All 
wires and electrical leads were thermally anchored to 
the top of the shield as well as the top-plate of the 
pressure-vessel. These modifications reduced the tem- 
perature difference between the top and the bottom of 
the resonator to less than 0.5 mK at all temperatures. 



Further discussion of temperature gradients appears in 
Sec. 5.3. 

The lowest temperature at which the gas-constant ap- 
paratus could be used for primary thermometry was 
determined by the 0-ring that was used to seal the body 
of the pressure vessel to its lid. When the temperature 
was too low, this seal leaked and either the methanol 
flowed from the bath into the pressure vessel or the gas 
from the pressure vessel flowed into the bath. The 
Viton^ 0-rings that had been used were replaced with 
silicone rubber 0-rings since they have a lower ultimate 
working temperature than their Viton equivalents. Even 
with this change, the 0-ring leaked at the lowest temper- 
ature of study, restricting the maximum pressure to 
360 kPa on the 217 K isotherm. The Viton 0-rings that 
sealed the transducers in their ports within the resonator 
were not replaced. Usually, the pressure difference 
across these 0-rings was only 1 kPa or less, with the 
higher pressure inside the resonator. If these 0-rings 
were leaking significantly during the intervals that the 
valve atop the resonator was closed, the resonance fre- 
quencies would have shown a systematic time depen- 
dence that was not detected. 

2.2 Dimensional Measurements 

The fabrication and characterization of the gas-con- 
stant resonator were described in detail in [4]. As dis- 
cussed in 1988 in Sec. 3.7 of Ref. [4], it was not possible 
at that time to obtain a satisfactory interpretation of the 
frequency- splittings of the non-radial (1,3) and (1,8) 
acoustic modes in terms of the geometry of the assem- 
bled resonator. (For the acoustic modes, the notation 
(/,n) is used to identify the order of the spherical Bessel 
function (/ = 0, 1, 2, . . .) and n denotes the number of 
nodes in the radial component of the acoustic velocity.) 
The problem persisted after measurements of the fre- 
quencies of microwave resonances in the same cavity 
[6,7] . We now reconsider this problem. 

The dimensions of each hemisphere have been mea- 
sured three times. The first series of measurements was 
made by the NBS shop that fabricated the hemispheres 
using their coordinate measuring machine. These mea- 
surements were completed before the resonator was as- 
sembled in 1985 and the resulting data are denoted 
R1985 in Table 3. In Table 3, three dimensions are given 
for each hemisphere: (1 ) as in Fig. 2 of Ref. [4], R^ is the 
exterior radius of the cylindrical surface near the equa- 
tor, (2) /?s is the average interior radius of the nearly 

' Certain commercial equipment, instruments, or materials are identi- 
fied in this paper to foster understanding. Such identification does not 
imply recommendation or endorsement by the National Institute of 
Standards and Technology, nor does it imply that the materials or 
equipment identified are necessarily the best available for the purpose. 
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Table 3. Dimensions of hemispheres (in mm) at 20. 1 "C 



Lower hemisphere 



Upper hemisphere 



Measurement series 


Re 


Rs 


hs 


Re 


Rs 


h 


R1985 


107.960 


88.935 




107.958 


88.918 




P1997 


107.967 


88.9133 


0.0585 


107.968 


88.8901 


0.0554 


F1997 




88.915 


0.056 




88.896 


0.041 



spherical surface, and (3) K is the length of the cylindri- 
cal extension of each hemisphere beyond the plane that 
would terminate a perfect hemisphere. 

After the resonator was disassembled in 1997, the 
second series of measurements, denoted PI 997, were 
made by the Precision Measurement Division of NIST. 
These measurements used their coordinate measuring 
machine, which was the most accurate machine avail- 
able to us. We were surprised by the comparatively large 
differences between the R1985 and P1997 results; 
therefore, we had a third series of measurements de- 
noted F1997 made by the Fabrication Technology Divi- 
sion of NIST using a third coordinate measuring ma- 
chine. (The machine used in 1985 was no longer 
available.) 

The F1997 measurements are consistent with the 
nominally more accurate PI 997 measurements. As indi- 
cated in Table 3, the average radius of each hemisphere 
7?s was approximately 0.02 mm smaller than reported in 
Ref. [4] and the average radius of the cylindrical boss R^ 
on each hemisphere was approximately 0.01 mm larger 
than reported in Ref. [4]. These differences are larger 
than the uncertainty of 0.005 mm claimed for the Rl 985 
data and we have no explanation for the differences. 
Fortunately, these dimensions were not critical for either 
the present study of ITS-90 or for the redetermination of 
the universal gas constant R. The cylindrical bosses 
were used to align the hemispheres during the assembly 
of the resonator. The R1985 and P1997 wlues for R^ 
differ; however, both series of measurements show that 
7?c of both hemispheres was nearly identical. This was 
confirmed by simple observations made with the hemi- 
spheres in contact with each other. 

The P1997 measurements did confirm that essentially 
all of the cylindrical extension of the hemispheres had 
been removed prior to the assembly of the resonator. 
This had been in doubt after the resonator was first 
assembled [4]. 

The disassembly of the resonator in 1997 revealed 
another surprise. When the transducers were inserted 
into their ports in the upper hemisphere of the resonator, 
the diaphragms of the microphones were not aligned 
with the interior surface of the hemisphere as intended. 
Instead, the transducers were recessed approximately 
0.9 mm in their ports. Thus, when the resonator was 



assembled, the interior of the upper hemisphere had two 
coin-shaped volumes 9.49 mm in diameter and 0.9 mm 
deep. Measured by total volume, these coin-shaped cav- 
ities were smaller departures from the intended spheri- 
cal figure than either the cylindrical extensions to the 
equators of the hemispheres or the difference between 
the radii of the hemispheres. The upper hemisphere and 
both transducer assemblies were machined in accor- 
dance with their drawings. We surmise that a design 
error was made, one surface of a 0.89 mm high locating 
step was used as the reference for designing the trans- 
ducer ports and, by accident, the other surface of the 
step was used as the reference for designing the trans- 
ducer assemblies. This error in design is similar to the 
error made during the ground tests of the mirror for the 
Hubble Space Telescope. 

The unintentional and unknown presence of the coin- 
shaped volumes did not compromise the redetermina- 
tion of R . For that work, the resonator's volume was 
determined by weighing the mercury required to fill it. 
During the weighings, plugs were substituted for the 
transducers and the coin shaped volumes were present 
and accounted for. See Ref. [4] for details. 

In summary, the cavity differed from a spherical fig- 
ure in three respects: (1) the two hemispheres had dif- 
ferent radii, (2) both hemispheres had cylindrical exten- 
sions, and (3) there were two coin-shaped recesses in the 
upper hemisphere. These departures from sphericity 
partially removed the degeneracy of the microwave 
modes [6] and of the non-radial acoustic modes [4] . 

Mehl [16], [17] expanded the shape of a cavity in 
spherical harmonics 



■■a[l - e^ '^c,,nYim(0,4>)], 

1=0 m=-l 



(8) 



and he provided explicit formulae for the splitting of the 
non-radial acoustic modes and for the second perturba- 
tions to the radial modes in terms of the expansion 
coefficients c;„ . Comparable results for the splitting of 
microwave triplets have also been derived [6]. The split- 
tings of the / = 1 microwave and acoustic modes depend 
on 6C20 only, and by far the largest contribution to 6C20 
comes from the cylindrical extensions. Here we com- 
pare the results of the present dimensional measure- 
ments with the splittings measured near 20 °C. 
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Source 

PI 997 dimensions 
F1997 dimensions 
(1,3) acoustic mode 
(1,8) acoustic mode 
TMl 1 microwave mode 
TM12 microwave mode 
TMl 3 microwave mode 



10^ec2, 

6.35 
5.40 
5.48 
5.63 
5.61 
5.56 
5.51 



Evidently, the observed splittings are more nearly con- 
sistent with the F1997 dimensions than the P1997 di- 
mensions. The agreement of these subtle features of the 
resonator with theory is encouraging. 



3. Microwave Measurements 

Microwave measurements were used to determine the 
volumetric thermal expansion of the spherical cavity in 
the temperature interval 213 K to 303 K. Measurements 
were carried out in the vicinity of the TMl 1 , TMl 2, and 
TMl 3 microwave triplets. Theory shows that the aver- 
age frequency of these triplets was insensitive to geo- 
metric imperfections that leave the internal volume un- 
changed [5]. In addition, the accuracy of microwave 
measurements was demonstrated by directly comparing 
microwave measurements with mercury dilatometry [6]. 

3.1 Apparatus and Data Acquisition 

For the measurements of /„, the acoustic transducer 
assemblies were removed from their ports and replaced 
with plugs from which the microwave coupling probes 
extended, exactly as described in [7]. These plugs had 
the same external dimensions as the acoustic transducer 
assemblies; thus, the ends of the plugs were recessed 
0.9 mm from the interior surface of the resonator. For 
most of the present measurements, the coupling probes 
protruded 4.07 mm beyond the ends of the plugs. After 
most of the data were acquired, a few measurements 
were made with the probes shortened to 2.31 mm and a 
few measurements were made with the probes extended 
to as much as 20 mm. The effects of changing the length 
of the probes are discussed below. 

During the microwave measurements, the interior of 
the resonator was evacuated while the surrounding pres- 
sure vessel was filled with helium to a pressure near 
10 kPa to facilitate thermal equilibration. The valve atop 
the resonator was left open. 

The temperature of the resonator was determined 
from two capsule-type platinum resistance thermome- 
ters that were embedded in metal blocks fastened to the 
top (thermometer #LN303) and bottom (thermometer 



#LN1 888002) of the resonator. The temperature as- 
signed to the resonator was always the average of the 
two thermometers. The thermometer's resistances were 
recorded after each frequency measurement so that each 
mode had its own unique determination of the average 
temperature. Immediately after completion of the mi- 
crowave measurements, the resistances of the ther- 
mometers were checked in a triple point of water cell. 
The values of R (r„) and the coefficients given in Table 
6 were used to calculate the temperatures for the mi- 
crowave data listed in Table Al. Although the intervals 
between complete recalibrations of these thermometers 
were quite long, the thermometers have a long history of 
good stability (see Table 5). We conservatively estimate 
that the uncertainty of the resonator's temperature dur- 
ing the microwave measurements was no more than 
2 mK relative to lTS-90 and this uncertainty propagates 
into relative standard uncertainties of [a (T)/a (T^)f and 
(T—Tgo) that are less than 0.1 X 10"^ because (1/a) 
da/dJis only about 16 X 10"'^ K"'. 

The microwaves were generated and detected by a 
Hewlett-Packard Model 8753B network analyzer which 
was connected to the resonator through a Hewlett-Pack- 
ard Model 85047A S-Parameter Test Set. The analyzer 
was configured to measure S12, which is defined (for 
properly terminated lines) as the complex ratio of 
voltage transmitted through the resonator to the voltage 
incident on the resonator. These instruments permitted 
measurements up to 6 GHz. We made detailed measure- 
ments in the vicinity of the nearly degenerate TMll, 
TM12, and TM13 triplets which occurred at 1.47 GHz, 
3.28 GHz, and 5.00 GHz, respectively. The triply degen- 
erate TE multiplets could not be detected with the 
straight probes that we used, and we decided to avoid the 
complications of dealing with multiplets with more than 
three components. In Ref. [6] it was shown that the 
thermal expansion determined from such highly degen- 
erate multiplets is consistent with that obtained from the 
TMl 1 and TMl 2 triplets. 

The frequency of the oscillator in the network ana- 
lyzer was derived from the highly stable quartz refer- 
ence oscillator installed in the audio frequency synthe- 
sizer used for the acoustic measurements. The frequency 
stability of this quartz oscillator was periodically con- 
firmed by comparison with primary standards. 

The network analyzer was used to scan 101 frequen- 
cies spanning each multiplet under study. The scan 
widths were 1.2 MHz, 1.6 MHz, and 2.0 MHz for the 
TMll, TMl 2, and TMl 3 triplets, respectively Typi- 
cally 20 scans with an IF bandwidth of 10 Hz were 
averaged. For each frequency, the averaged values of the 
real (u) and the imaginary (v) parts of the signal trans- 
mitted through the resonator were down-loaded to a 
computer for fitting. The values of u and v at all of the 
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101 frequencies were fit to a sum of two Lorentzian 
functions of the frequency, 



U + IV- 



1 



ifA„ 



r - Fl) 



jr + B + C(f-f,). 



(9) 



Here, A,B, and C are complex constants, and F^ =fm + 
igm are the complex, nearly-degenerate resonance fre- 
quencies of the triplet under study. In Eq. (9), the 
parameters B and C account for possible crosstalk and 
for the effects of the "tails" of the modes other than the 
one under study. Although the multiplets we studied 
were expected to be nearly degenerate triplets, the data 
could be represented by Eq. (9) summing over just two 
complex resonance frequencies (see Fig. 3). Typically 
the standard deviation of the fit, expressed as a fraction 
of the maximum amplitude measured, was 0.00016, 
0.00035, and 0.00075 for the TMll, TM12, and TM13 
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Fig. 3. Top: Measurements of the in phase and quadrature signals 
detected as the microwave generator was swept through the TMll 
"triplet." Bottom: Deviations of the detected signals from a fit of Eq. 
(9) to the data in the upper panel. In this particular case, the fit 
determined 10 parameters: two complex resonance frequencies, two 
complex amplitudes, and a complex additive constant. 



triplets, respectively. Thus, the present signal-to-noise 
ratio was at least a factor of 4 larger than that achieved 
in Ref. [6]. The repeatability of the fitted frequencies 
and half- widths were always less than 100 Hz (or equiv- 
alently, (0.02 to 0.07) X 10~^ of the resonance frequen- 
cies). The frequencies and half-widths resulting from 
the fits to the microwave data are in Table Al in the 
Appendix. 

For the TMl 1 and TM13 triplets, the deviations from 
the fitted functions were random. This was confirmed 
by further averaging. The deviations were reduced by a 
factor of 3 and remained random. For the TM 1 1 triplet, 
B and C in Eq. (9) were negligible. For the TM 1 2 triplet, 
the deviations from the fitted functions were systematic 
and the parameter B in Eq. (9) accounted for nearly one 
quarter of the peak signal. These observations are a 
consequence of the proximity of the TM4 1 multiplet to 
the TMl 2 triplet and the comparatively efficient cou- 
pling of the TM41 multiplet to the probes. (The peak 
amplitude in the vicinity of the TM41 multiplet was 30 
times the peak amplitude in the vicinity of the TMl 2.) 
Thus, it is not surprising that we were unable to improve 
the quality of the fit to the TMl 2 triplet by adding a 
third resonance frequency within the scanned range 
with a half-width comparable to the other two. 

3.2 Interpretation of Microwave Frequency 
Measurements 

3.2.1 Partial Splitting of the Triplets 

In Refs. [6,7], as in this work, the TMll and TM12 
triplets were fitted by just two resonance frequencies. In 
Ref. [6], it was argued that the third component of the 
expected triplet was not detected because the deviations 
of the cavity from a perfect sphere were primarily 
axisymmetric. Furthermore, it was argued from acous- 
tic phase measurements that the lower frequency com- 
ponents of the TMl 1 and TMl 2 triplets were unresolved 
doublets and the upper components were singlets. We 
obtained further evidence supporting this conclusion by 
soldering a 20 mm long curved wire to one of the 
probes. This extension produced a partial splitting of all 
three components of the triplets. A multi-parameter fit 
to the TMl 3 triplet resulted in a splitting of 250 kHz 
between the lowest two components and a splitting of 
602 kHz between the highest two components. The lat- 
ter is close to the 660 kHz splitting between the doublet 
and the singlet that was observed with short probes. The 
fractional splittings [8/7(/m) in Eq. (10)] range from 
111 X lO"*^ to 222 X 10-*^ (see Fig. 4). 
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Fig. 4. The temperature dependence of the frequency differences 
between the two resolved components of three microwave triplets, 
expressed as a fraction of the average frequency of each triplet. The 
open symbols denote data taken with shortened probes. 

"flifm) — (/singlet ~ /ioubletVCj /singlet + j/ioublet)- (10) 

Following Ref. [6] we note that if the splitting resulted 
from an axisymmetric deformation of the spherical cav- 
ity's radius a such that the radial coordinate r was given 
by 

r = a[l-e2cmYio{e)l (11) 

1=0 

the splittings of the 1=1 TM triplets would be 






ec2o 



(Vn,)' 



(12) 



(Vm is a microwave eigenvalue.) We have ignored any 
possible deformations with l>2 and used the measured 
splittings together with Eq. (12) to determine ec2o- For 
the three microwave triplets we found ec2o = — 560 X 
10~''. Remarkably, the values of ec2o differ from their 
mean by less than 1 % at each temperature. (See Fig. 5.) 
The negative values of ec2o imply that the polar "radius" 
of the resonator fractionally exceeds the equatorial 
"radius" by approximately 560 X 10~''. The tempera- 
ture dependence of ec2o indicates that the thermal expan- 
sion of the resonator is not isotropic. As the temperature 
is increased from 213 K to 303 K the average radius of 
the cavity increases fractionally by approximately 
1.4 X 1 0~^. However, the fractional increase of the polar 
radius is 5 X 10~* less than that of the equatorial radius. 
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Fig. 5. The temperature dependence of the deformation parameter 
—€C2o calculated by substituting the data of Fig. 4 into Eq. (12). Note 
that the zero of the ordinate is suppressed and that the data for the 
three modes are nearly consistent. This is evidence that the microwave 
splittings are mostly the result of a simple shape imperfection. The 
open symbols denote data taken with shortened probes. 

The open symbols in Fig. 5 indicate the values of ec2o 
obtained when the length of the probes was reduced 
from 4.07 mm to 2.31 mm. It is apparent that the probes 
had minor but detectable influence on the observed 
splitting of the TMl 1 mode. 

3.2.2 Widths of Microwave Resonances 

The penetration of the microwave field into the wall 
of the resonator results in a contribution to the half- 
widths of the resonances g^ as well as an equal reduc- 
tion in the frequencies of the resonances. The magnitude 
of these effects for the TMlm modes considered here is: 



A/m +igm , . , ,-A 

/n, ^ '2a 



2_ 



(13) 



In Eq. (13), 8^ is the electromagnetic field penetration 
length given by 



S„ = Viif^ficr, 



(14) 



in which a is the conductivity of the stainless-steel shell 
and fi is the magnetic permeability. 

In order to calculate g^, we assumed that the mi- 
crowave permeability of the cavity's wall was exactly 
fio, the permeability of free space, and that the electrical 
conductivity of the cavity's wall at microwave frequen- 
cies was identical with the dc conductivity. The latter 
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was estimated from the data of Clark et al. [18] who 
measured the resistivity of samples of 316 stainless- 
steels from four different manufacturers at five tempera- 
tures: 273 K, 192.4 K, 75.75 K, 19.65 K, and 4.0 K. For 
the three highest temperatures, the resistivity of each 
sample is very nearly a linear function of temperature. 
We represented the resistivity data in the range 2 1 3 K to 
303 K by the following linear function of temperature T: 

p/(li • m) = 7.522 X 10""[1 + 1.26 X 10"'(r- To)], 

(15) 

where To = 273. 15 K. At 273 K resistivities of four sam- 
ples ranged ±4 % about the value returned by Eq. (15). 
The corresponding range in calculated electromagnetic 
penetration lengths is ±2 %. Using Eqs. (13) to (15), 
we calculated g„(calc.)//; = (87.2 ± 1.7) X 10"*^ for the 
TMl 1 triplet and g„(calc.)//„ = (29.8 ± 0.6) X lO"*^ for 
the TM13 triplet. The uncertainty in gn,(calc.) propa- 
gates into an uncertainty of 0.39 X 10 ^ in [a{T)l 
a(T'„)]^ at the lowest temperature of this study and 
0.18 X 10 "^ to the uncertainty in [a{T^la{T^)f. 

Figure 6 compares the experimental values of g^ to 
the calculated values of g^. The plot shows the scaled 
excess half-width Ag^ defined by 



Agm//m = [gm(eXpt.) - gn,(calc.)]//n, 



(16) 
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Fig. 6. Temperatures dependence of the fractional excess half-widths 
of the single (s) and doublet (d) components of three microwave 
triplets [Ag„ = g„(meas.) — g„(calc.)]. The open circles represent the 
data taken with the shortened probes. The dashed curves show the 
same data recalculated with the assumption that the electromagnetic 
field penetration length is increased by 10 % over that calculated from 
the low-frequency resistivity. 



The values of Agjfin are positive, nearly temperature-in- 
dependent, and differ for each component of each 
triplet. (In Fig. 6, the higher frequency components are 
labeled s and the lower frequency components are la- 
beled d.) If Agjfin had a temperature dependence as 
strong as that of gn,(expt.)//n, it would have been de- 
tected. 

Figure 6 shows that the effect of reducing the length 
of the coupling probes from 4.07 mm to 2.31 mm was 
detectable, but small. [Compare the solid circles (long 
probes) with the open circles (short probes).] Thus, the 
probes cannot explain most of Agm(expt.)//in. 

We noted that Agjfin varied approximately as (/m) "^■ 
This suggested that the phenomenon responsible for the 
excess half- widths was an evanescent wave and led us to 
reconsider the half-width data obtained with the short- 
ened probes. We made the ad hoc assumption that the 
electromagnetic field penetration lengths were 10 % 
longer than those calculated above. This would be the 
case if the resistivity of our shell were 20 % larger than 
that measured for comparable alloys by Clark et al. [18] 
or if the magnetic permeability of our shell were 20 % 
larger than ixo- (The relative magnetic permeability of 
type 316 stainless steel is reported as 1.008. Cold work- 
ing increases the permeability of type 316 slightly and 
greatly increases the permeability of similar alloys with 
slightly lower nickel content [19].) Using the larger pen- 
etration lengths, we found that Agjfin for all of the 
microwave components fell in a narrow range spanning 
zero (see the dashed curves on Fig. 6). With this same 
ad hoc assumption, the volumes determined from the 
three triplets span a range of only 4 X 10" '', fractionally. 
This contrasts with the range of 18 X 10"^ deduced 
from the smaller penetration lengths. The consequences 
of this alternative analysis on T/T^ were calculated and 
the difference between the two analyses contributed to 
the uncertainty of T/T^ that appears in Row 3 of Table 
2. 



3.3 Volumetric Thermal Expansion of the 
Resonator 

To determine the thermal expansion of the resonator, 
we computed a "corrected" weighted average radius 
(flm) for each triplet at each temperature. Because the 
lower frequency component of each triplet is an unre- 
solved doublet, we defined the average by 



{a„) = 



'2'n[-^fdoub\et + 3./singlet + gm(calc.)] 
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in which gj^calc.) is the value of the half-width com- 
puted from Eqs. (13), (14) and (15). With the probes at 
their "normal" 4.07 mm length, the values of (oni) for the 
TMl 1 triplet were always about 10 X 10"'' larger than 
the values of («„) for the other triplets. (When the 
probes were shortened, this discrepancy was greatly re- 
duced.) 

The microwave frequencies were measured during the 
1989 runs and spanned the temperature range 210.74 K 
to 302.93 K. Because the microwave frequencies had not 
been measured near the 217 K isotherms used for the 
1992 acoustic measurements, an interpolation function 
for the thermal expansion was needed. We fitted a poly- 
nomial function of (T^—T) to the three triplets of 
{aJ(T)). The results can be represented by 



10' 



{am{T)) V 



1 



+ 0.0170417" + 0.0000411' 



= -31.3144371 



(18) 



with 1= (T^—T)/K. All the parameters in Eq. (18) had 
a significance greater than 0.999 based on the F-test and 
the Student ^-distribution. The fractional standard devia- 
tion of the fit for [{a(T))/{a(T^))f was 0.22 X 10"". 
The deviations are shown in Fig 7. Clearly, the measure- 
ments of the microwave frequencies define a thermal 
expansion function with extraordinarily high precision. 
The curved line in Fig. 7 shows the effect of replacing 
the calculated halfwidths gni(calc.) in Eq. (17) with 
gni(expt.), the experimental values. The effect of this 
change is slight. The choice of the functional form for 
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Fig. 7. Fractional deviations of the microwave data from the temper- 
ature-dependent average internal radius (u„) of the resonator. The 
average radius for each TMlm microwave triplet is defined by Eq. (17) 
with the parameters of Eq. (18). The solid curve shows the small effect 
of replacing the theoretical values of the half-widths in Eq. (17) with 
their experimental values. 



Eq. (18) is not critical. Indeed, if linear interpolation 
between adjacent microwave isotherms 210 K and 
224 K had been used to obtain [{aiT))/{aiT^))f on the 
acoustic isotherm 217 K, the result would have differed 
by only 0.48 X 10"" from Eq. (18), and this is the worst 
case. 

We considered several alternatives to defining the 
radius by Eq. (17) and fitting {a^{T)) by Eq. (18). In 
reviewing the alternatives, recall that increasing [a{T)l 
a{T^)f by 1 X 10"" increases TIT^ by 1 X 10"" or, 
equivalently, 0.30 mK at T^. If we had used the mea- 
sured half-width of each component of the multiplet 
instead of the calculated half -width in Eq. (17), [a{T)l 
a{T„)f would have changed by 0.24 X 10"" in the 
worst case. If we had not deduced that the lower compo- 
nent of each triplet was an unresolved doublet and used 
the definition 

{P'm) = 7^ rr~7 I r~7 j 7 1 vT ' ( ^ ^) 

^'"■[2 /doublet + 2/singlet + gra(,CalC.JJ 

the change in [a (T)/a (T^)]^ would have been less than 
0.06 X 10"". If we fit the data for each multiplet sepa- 
rately, the range of the values of [a (T)/a (T^)f would 
have been 0.4 X 10"". The small effects of these 
changes is reassuring. 



3.4 Comparison with Earlier Thermal Expansion 
Measurements 

The present measurements of the volumetric thermal 
expansion are compared to earlier microwave measure- 
ments and to mercury dilatometry in Table 4. Using 
the weighted average of the TMlm modes we find 
10"[yo(rg)/yo(rw) - l] = 1419.0 ± 0.7. The discrepancy 
with the previous measurements is displayed in Table 4 
and is equivalent to 0.25 X 10 " in Vo(T^yVo(T„), when 
the ratio from Ref. [6] is recalculated using the values 
of gjn obtained from Eqs. (13) to (15). (In Ref. [6] the 
expression used for the resistivity differed slightly from 
Eq. (15).) The two sets of microwave measurements 

Table 4. Volumetric expansion of the cavity between T„ and Tg 



Mode 


10''[ViTg)/ViT^)- 


-1] 


Reference 


TMll 


1418.7° 
1418.7 






[6] 
This work 


TM12 


1418.9° 
1419.3 






[6] 
This worl£ 


TM13 


1419.2 






This work 


Mercury dilatometry 


1416.6± 


1.5 




[6] 



'Recalculated using g„ from Eqs. (13), (14), and (15). 
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agree to well within estimated uncertainties. The frac- 
tional difference between this determination and the 
mercury dilatometry is 2.45 X 10 ''; this is within 1.5 
combined standard uncertainties. 



4. Acoustic Measurements 

With the minor exceptions mentioned below, the 
complex resonance frequencies F^=f^+ ig^ of the radi- 
ally symmetric, non-degenerate acoustic modes were 
determined with the same transducers, instruments, and 
procedures that were used with this resonator in the past 
[3,4,7]. For brevity, we describe the exceptions and the 
changes while omitting repetition of details that are 
unchanged. 

For argon, the frequencies and half-widths of the 
lowest five radial modes [conventionally designated 
(0,2), (0,3), . . . (0,6)] were measured. As discussed 
below, the (0,6) mode had a small, anomalous pressure 
dependence that may have resulted from a near coinci- 
dence of the frequency of that mode with a mode of the 
spherical shell. To explore the effect of this, a separate 
analysis was conducted omitting all the argon data for 
the (0,6) mode. 

For xenon, the frequencies of lowest seven acoustic 
modes were measured; however, only the lowest five 
were used in the final analysis. We knew (see Fig. 8 of 
Ref. [4]) that the resonance frequencies of the (0,7) and 
(0,8) modes partially overlap those of nearby, highly 
degenerate, non-radial modes; however, we expected 
that the effects of the overlap could be canceled to a high 
degree when computing speed-of-sound ratios on a 
mode by mode basis, as implied by Eq. (7). However, 
the uncertainties of the xenon data were dominated by 
unanticipated impurity effects. This, and the fact that 
the mode-by-mode analysis proved more cumbersome 
than anticipated led us to discard the data for the (0,7) 
and (0,8) modes. 

In Ref. [4], arcing within the detector transducer was 
reported. On occasion, the problem of arcing within the 
detector transducer reappeared. It was permanently 
cured by reducing the dc bias voltage on the detector 
transducer from 200 V to 150 V. The reduction of the 
bias voltage reduced the signal-to-noise ratio approxi- 
mately 25 %. This was more than offset by altering the 
protocol for measuring the signal produced by the de- 
tector transducer. In our earlier work, an analog lock-in 
amplifier had been used to measure the detected acous- 
tic signal and the amplifier's post-detection filter was 
used for averaging. Then, one had to wait eight filter 
time constants for the output to settle fractionally to 
within 3 X lO""* of its final value before recording the 
output, and the recorded output benefitted from one-fil- 



ter-time-constant of averaging. In this work, we used a 
digital lock-in amplifier. The post-detection time con- 
stant of this amplifier was set to 0.3 s. After each incre- 
ment of the acoustic frequency, we waited for the output 
of the lock-in amplifier to settle to within 10""* of its 
final value. (The settling occurred with the time con- 
stant Ts = 1/(217^2)) Then, the output of the lock-in am- 
plifier was measured at intervals of Ts and digitally aver- 
aged. In comparison with our previous work with an 
analog lock-in amplifier, the measurement time was un- 
changed and the signal-to-noise ratio was increased up 
to a factor of 3 as the pressure was reduced below 
100 kPa. (In previous work [4] below 100 kPa, Ts < 1 s, 
the post-detection filter had been set to 3 s, and the 
signal-to-noise ratio decreased as p~^.) 

For both argon and xenon, each complex acoustic 
resonance frequency Fa was determined by fitting the 
detector's response to a single Lorentzian function of 
the frequency (Eq. (9) with the summation index m 
taking on the value 1, only). Equation (9) was fitted to 
the data for each mode twice, once with the constraint 
C = 0, and once without the constraint. The constrained 
fit was used for further analysis, except when relaxing 
the constraint reduced the standard deviation of the fit 
by at least 30 %. 

For weighting fits of the temperature and pressure 
dependencies of /a, it is convenient to have an approxi- 
mate expression for the standard deviation of a single 
measurement of /a. An expression for argon was devel- 
oped during the redetermination of R and of Tg. 



(r(f,) = 10-ya[l + (100 kFa/pf{6 kUz/f^f 



(20) 



At pressures above 100 kPa the signal-to-noise ratio was 
sufficiently high that the imprecision of a measurement 
was dominated by small, uncontrolled phase shifts in 
the measurement system. The loss of precision at low 
pressures is a consequence of the signal declining asp^'^ 
and the resonance half-widths increasing as p~"^. When 
the resonator was filled with xenon the transducers' 
characteristics were essentially unchanged. For xenon, 
Eq. (20) is a reasonable estimate for (r(f^), provided that 
the characteristic pressure in that equation is replaced 
with 55 kPa, in accordance with expectations based on 
xenon's thermophysical properties. In this work, we 
continued to use Eq. (20) with the characteristic pres- 
sures 100 kPa and 55 kPa to weight the data, even 
though Eq. (20) over-estimated the uncertainty of the 
lowest-pressure data. 

The systematic errors arising from the reference os- 
cillator in the frequency synthesizer, non-linear effects, 
and the instrumentation for frequency measurement 
were found to be negligible. 
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Prior to the commencement of the acoustic measure- 
ments, the resonator was evacuated and "baked" at 
50 °C for 4 d. Simultaneously, the gas manifold was 
baked at or above 100 °C. At the conclusion of the 
baking procedure the residual pressure indicated by an 
ionization gauge was approximately 13 [xPa. After the 
acoustic measurements at r„, T^, and 217 K, the res- 
onator was again baked, this time at 60 °C for 24 h. 
After the conclusion of the argon measurements and 
prior to the xenon measurements at T^, the gas manifold 
was baked again at 1 00 °C for 24 h to remove any 
residual traces of argon from the pipework. The res- 
onator was baked again before the final xenon measure- 
ments at Tg. 

At T^, the acoustic measurements were made at 13 
pressures between 25 kPa and 500 kPa in three separate 
fillings of the resonator. On the other isotherms, the 
measurements were made at pressures which corre- 
sponded to approximately the same densities. (Only at 
217 K, where the 0-ring sealing the pressure-vessel 
leaked, was it necessary to limit the maximum working 
pressure.) In this way, we ensured that the number of 
parameters required to fit the virial expansion to the 
acoustic data was the same on every isotherm studied 
and we avoided the possibility of biasing TIT^ by using 
different orders of fit on various isotherms. In xenon, 
measurements were conducted at pressures in the range 
30 kPa to 300 kPa. 

For both argon and xenon, the frequencies and half- 
widths of the radially symmetric modes were deter- 
mined twice at each temperature and pressure, first in 
order of ascending mode index, and then in the reverse 
order. The resistances indicated by the three thermome- 
ters were recorded after each frequency measurement, 
ensuring that each mode had a unique determination of 
the mean temperature of the cavity. 



5. Thermometry 



5.1 Temperature of the Bath 

A calibrated thermometer was used to measure the 
temperature of a point within the bath. This facilitated 
adjusting the bath's temperature to within a few mK of 
the resonator's temperature. Because the bath's temper- 
ature was so close to the resonator's temperature, the 
latter usually drifted less than 0.5 mK during the 45 min 
required to acquire the acoustic data at each temperature 
and pressure. 



5.2 Temperature of the Resonator 

The temperature of the gas within the cavity was 
inferred from three capsule-type standard platinum 
resistance thermometers (SPRTs). Two of these ther- 
mometers (LNl 888002 and LN303) had been used dur- 
ing the redetermination of R and of Tg. These were 
mounted in metal blocks attached to the bosses at the top 
(LN303) and bottom (LNl 888002) of the resonator. 
The third SPRT (Chino serial number RS18A-5) was 
purchased after the microwave measurements were 
completed and installed near the equator of the resonator 
prior to the acoustic measurements of 1989. This ther- 
mometer was surrounded by two metal sleeves that were 
threaded into the resonator. The inner sleeve was OFHC 
copper and the outer was aluminum. The resistances of 
each thermometer were measured before and after each 
measurement of /s and the temperature associated with 
that measurement was the average of the temperatures 
of the three thermometers. 

5.3 Thermometer CaUbration, History and 
StabUity 

We now describe the other factors which lead to the 
uncertainty estimates in Table 2 under the heading 
"Thermometry". The resistance bridge, its standard re- 
sistors, and the SPRTs together function as a transfer 
and interpolation standard between the triple point cells 
and the resonator. Therefore, the primary concern is the 
long term stability of the thermometers. This was evalu- 
ated by periodically checking the thermometers in triple 
point cells. Table 5 lists the quantity R{Ti, i — > 0) which 
is the resistance measured with the 30 Hz ac bridge 
extrapolated to zero current. Here, T", represents the tem- 
perature of the triple point t where the subscript "t" 
stands for "w", "m", "g", or "a" in the notation r„, T^, 
Tg, or 71 which denote the temperatures of the triple 
points of water, mercury, gallium, or argon, respectively. 
(In contrast, Refs. [4] and [7] use the symbol T", to 
represent the temperature of the triple point of water.) 
The average change in R{Ti, i — > 0) between our cali- 
brations is 10 (xfi. We used this value as the estimated 
standard uncertainty in (T—Tgo)/T resulting from our 
imperfect thermometer calibrations and listed the results 
in Row 8 of Table 2. 

The calibrations at T„, T^, and Tg just before (May 
1992) and just after (August 1992) the acoustic mea- 
surements, together with subsequent (October 1992) 
calibrations at r„ and 71 were used to generate two sets 
of coefficients for the defining equations on ITS-90. The 
two sets of coefficients were used to calculate two tem- 
peratures (Tbeforc and TkfEr) for each isotherm. The tem- 
perature differences (Tbrforc— Tlfter) are one measure of 
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Table 5. Record of calibration of thermometers 



Table 6. Summary of thermometer characteristics 



Date 



07/04/92 
28/05/92 
27/08/92 
23/10/92 

07/04/92 
28/05/92 
27/08/92 
23/10/92 

28/05/92 
27/08/92 
23/10/92 



Thermometer 



t = Water 



LN1888002 
LN1888002 
LN1888002 
LN1888002 

LN303 
LN303 
LN303 
LN303 

RS18A-5 
RS18A-5 
RS18A-5 



t = Gallium 



t = Argon 



R(Tt,i ^oyn 



25.541502 
25.541497 
25.541518 
25.541527 

25.475522 
25.475521 
25.475542 
25.475532 

25.206136 
25.206148 
25.206148 



28/05/92 


LN1888002 


28.558768 


27/08/92 


LN1888002 


28.558760 


28/05/92 


LN303 


28.484839 


27/08/92 


LN303 


28.484831 


28/05/92 


RS18A-5 


28.182858 


27/08/92 


RS18A-5 


28.182863 


t = Mercury 


28/05/92 


LN1888002 


21.560995 


27/08/92 


LN1888002 


21.561004 


28/05/92 


LN303 


21.505550 


27/08/92 


LN303 


21.505572 


28/05/92 


RS18A-5 


21.279150 


27/08/92 


RS18A-5 


21.279155 



23/10/92 


LN1888002 


5.515220 


23/10/92 


LN303 


5.502255 


23/10/92 


RS18A-5 


5.449206 



the uncertainty of T resulting from imperfect calibra- 
tions. The scaled differences 10*' X (rbeforc~7"after)/rare 
listed in Row 9 of Table 2. These differences arise from 
a common source; thus, they are not random. Indeed, 
they are highly correlated. 

Table 6 lists values for the parameters a and b which 
occur in the definition of the International Temperature 
Scale of 1990. These values were determined from the 
calibrations before the acoustic measurements, given in 
Table 5, in conjunction with the triple point of argon 
measurements performed by the NIST Thermometry 



Thermometer 


LN1888002 


LN303 


RS18A-5 




T/K S T„ 


(May 1992) 




lO^a 
10*^6 


-8.545941 
2.636814 


-14.787059 
3.918078 


-40.359244 
6.610082 




Tw < T/K S 


Tg (May 1992) 




lO^a 


-1.245599 


-1.774805 


-4.376422 



T/K S T„ (August 1992) 



lO'a 
lO'^fc 



-8.299809 
4.242863 



-14.871504 
3.367348 



-40.211829 
7.572838 



Tw < T/K S T, (August 1992) 



lO'a 



-1.350379 



-1.879442 



-4.406030 



Group after the completion of the acoustic measure- 
ments. 

To determine temperatures from the experimental 
resistance ratios, W(Tgo), the specified deviation func- 
tions for the appropriate ranges were used in conjunc- 
tion with the parameters given in Table 6 to determine 
the reference function Wr(Tr)o)- The specified inverse 
function was then used to calculate approximate tem- 
peratures on Tr)o- The derivative of this function with 
respect to WrCTgo) was determined numerically and used 
to adjust the Tgo temperatures such that Wr(Tgo) calcu- 
lated from the deviation function and the reference 
function in conjunction with the assigned Tgo value 
agreed to better than 10~'. This ensured that the Tgo 
temperatures determined from Wr(Tr)o) could be calcu- 
lated with arbitrary precision for a given resistance ratio 
W(Tgo)- In this way, any error in the calculated temper- 
atures resulted entirely from errors in the calibrations, 
sub-range inconsistencies, and non-uniqueness of the 
scale. The first of these errors is estimated from re- 
peated calibrations at the various triple points, and the 
last two are inherent limitations of the scale for ther- 
mometers having equivalent calibrations. Above, we 
have specified the sub-ranges that we have used and 
have indicated them in Table 6. Thus the sub-range 
inconsistences do not contribute to the uncertainties of 
(T—Tr)o) if "Tgo" is understood to mean the specified 
sub-range. 

The non-uniqueness of ITS-90 does not contribute to 
the uncertainty of our values for either Tg/T^ or TJT„; 
however, it does contribute to the uncertainties of 
(T—Tr)o) on our isotherms at 217 K, 253 K, and 293 K. 
From Fig. 1 .5 of Ref. [20], we estimate the non-unique- 
ness of ITS-90 is 0.2 mK at 217 K and 253 K, provided 
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the thermometers are calibrated at r„, T^, and T^. Be- 
cause of the dearth of relevant published data, we made 
two naive observations to estimate the non-uniqueness 
contribution to the uncertainty of {T— Tgo) near 293 K. 
First, we noticed that the non-uniqueness varies approx- 
imately parabolically over the interval between closely 
spaced calibration points, and second, the interval 
(T^—T^) is approximately 0.77 X (r„— 7"^). Thus, we 
expect the contribution to be (0.77)^ times the 0.2 mK 
contribution at 253 K. The non-uniqueness contribu- 
tions to the uncertainty of (T—Tgo)/T appear in Row 1 1 
of Table 2. 

5.4 Calibration Techniques 

The resistances of the thermometers were measured 
with a four-wire, ac resistance bridge operated at 30 Hz. 
The bridge was designed by R. D. Cutkosky and built at 
NBS and has been designated NBS/CAPQ Microhm 
Meter 5. The bridge was used for all of the thermometer 
calibrations and all of the temperature measurements; 
thus, uncertainties associated with differences between 
bridges were not present [21]. This bridge is the same 
one which had been used for the re-determinations of R 
and Tg [4,7]. Normally, the bridge was operated with a 
measurement current of 1 mA and, with our 25 fl ther- 
mometers installed in the resonator, a typical standard 
deviation of reading was 3 |jufl. 

For calibration, each capsule thermometer was in- 
stalled in an extension probe similar to that described in 
Ref. [22]. Calibration resistance measurements were 
performed with currents of 1 mA and 2 mA so that we 
could extrapolate the results to zero current. The calibra- 
tion measurements were always performed in the order 
1 mA, 2 mA and 1 mA during a period of at least 10 
minutes after the SPRT had equilibrated in the triple- 
point cells. When unexpected drifts occurred, they were 
observed and the cause was eliminated prior to repeat- 
ing the measurement. In this way a standard deviation of 
a calibration measurement was 1 |jufl. When installed in 
the resonator, the self -heating of the thermometers was 
within 10 (xfl-mA~^ of that measured when the ther- 
mometers were installed in the calibration probes. 

The gallium point cell used for the calibrations was 
the same one that had been used in the redetermination 
of Tg. As mentioned in Ref. [7], B.W Mangum (then in 
the Temperature and Pressure Division of the NBS) 
compared this cell with one of a group which he main- 
tained as standards. The two cells were indistinguishable 
at the level of 50 |juK. For storage this cell was filled with 
argon; it was evacuated prior to use. 

The mercury point cell was manufactured and filled 
by us for this project following the guidelines of 
Furukawa et al [23]. The mercury used in the cell was 



a portion of the sample that had been used to determine 
the volume of the resonator in connection with the 
redetermination ofR. This mercury is traceable to A. H. 
Cook and came from the same NBS stock as described 
in Refs. [4] and [7]. All glass used in the construction 
of the cell was cleaned with HF prior to use. 

In March 1998, G. F. Strouse of NIST's Thermome- 
try Group made two direct comparisons of our mercury 
triple-point cell (Hg 88-1) with the Thermometry 
Group's laboratory standard mercury triple-point cell 
(Hg SS-1). The measurement system included an ASL 
Model F18 bridge operating at a frequency of 30 Hz 
with a 100 n Tinsley Model 5685 reference resistor and 
a 25.5 n SPRT. Corrections were made for the differ- 
ence in hydrostatic head effects due to the different 
immersion depths. After corrections, the triple-point 
temperature of the Hg 88-1 cell was 50 (xK lower than 
that of the Hg SS-1 cell. A standard uncertainty of 
0.10 mK was attributed to the value of T^ realized in the 
Hg SS-1 cell to account for the impurities and measure- 
ment uncertainties. 

The calibration at the 7"^ reported in Table 5 was 
made by G.F. Strouse of the Thermometry Group of 
NIST using the resistance bridge mentioned above. 

Our realization of the metal triple points was by 
means of the "double-melt" method. This method was 
selected in preference to the freezing method because it 
circumvents the problem of the massive undercool often 
observed when performing a "double-freeze." 

5.5 Temperature Gradients in the Resonator 

During the acoustic measurements a small vertical 
temperature gradient existed in the resonator. (See Fig. 
8.) The gradient was nearly symmetrical about the 
equator. In the worst case, when the bath was at 217 K, 
the top of the resonator was 0.5 mK warmer than the 
bottom of the resonator. The gradient may have been 
caused either by an undetected heat leak from the labo- 
ratory to the top of the resonator or by temperature 
stratification within the stirred bath surrounding the 
pressure vessel. Either possible cause is consistent with 
the observation that as the temperature was increased 
from 217 K through room temperature to Tg, the magni- 
tude of the gradient decreased, reaching zero at room 
temperature, and then increased with reversed sign. 

One can show that the acoustic resonance frequencies 
/a for the radially symmetric modes are determined by 
the volume average of the temperature distribution 
within the gas. Furthermore, if the deviation of the 
temperature from the mean is represented by 
8r(r,6),(/)) the frequency shifts are of the order of 
(8r^(r,6', (/))). Thus, only the asymmetry in the mea- 
sured vertical gradient contributes to the uncertainty of 
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T; this was probably less than 0.1 mK and therefore 
negligible. 

We had no means of detecting a horizontal compo- 
nent to the temperature gradient. If such a gradient were 
present, it probably was smaller than the vertical gradi- 
ent because the resonator, the cylindrical copper heat 
shield enclosing it, and the pressure vessel enclosing 
both were all axisymmetric (Fig. 2). Furthermore, the 
dimensions of the heat shield were such as to screen 
horizontal gradients more effectively than vertical gradi- 
ents. Somewhat arbitrarily, we assumed that any asym- 
metry of the horizontal temperature gradient was one 
half the temperature difference between either end of 
the resonator and its middle. Row 10 of Table 2 lists the 
uncertainty in (T—Tgo)/T corresponding to V2 the tem- 
perature difference between the middle of the resonator 
and the top of the resonator. 

As observed in Ref. [4], closing the isolation valve in 
the top of the resonator heated the thermometer at- 
tached to the top of the resonator 1 mK to 3 mK. Before 
measuring the acoustic frequencies, we waited until this 
thermal transient had decayed and the thermometers 
indicated temperatures that were identical with the ones 
before the valve was closed. 



6. Gas Samples 

The argon was withdrawn from the cylinder that con- 
tained the gas used in the re-determination of R and 7g. 
This cylinder had been purchased from Matheson Gas 
Products. The supplier's lot analysis provided the fol- 



lowing upper bounds for the mole fractions of impuri- 
ties: N2<3 X 10"^ 02< 1 X 10"'^; H20< 1 X 10'*^; 
and total hydrocarbons < 0.5 X 10"''. No further purifi- 
cation was attempted. The small effects of these impuri- 
ties on the speed of sound in argon can be found in Ref. 
[4]. 

The xenon was "research-grade" purchased from 
Matheson Gas Products. The supplier's lot analysis 
provided the following upper bounds for mole fraction 
of impurities: N2<2X10"^ O2 < 1 X 10"^ and 
Kr < 18 X 10"^. This sample was purified by exposing 
it for 91 h to a zirconium-aluminum alloy getter main- 
tained at 673 K. Under these conditions, the getter is 
very effective at removing gases such as CO, CO2, O2, 
N2, and H2 from the noble gases and it is moderately 
effective at removing hydrocarbons. After purification, 
the xenon was condensed into a small stainless steel 
cylinder for storage prior to being admitted into the 
resonator. This storage cylinder had been baked under 
high vacuum at a temperature above 400 K until the 
residual pressure was less than 10 (xPa (about 10~' torr). 
The change of the speed of sound in xenon upon addi- 
tion of an impurity with a mole fraction x « 1 can be 
calculated from the quantity (l/u^)(du^/dx). For the im- 
purities of interest, this quantity has the wlues: N2, 0.52; 
O2, 0.49; H2O, 0.43; CO2, 0.40; Ar, 0.70; and Kr, 0.36. 
The monatomic impurities that were not removed by the 
getter would not affect the ratios T/T„ unless their mole 
fractions changed (for example, by preferential adsorp- 
tion) when the temperature changed. (See Sec. 8.4 for a 
detailed discussion of possible contamination.) 
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7. The Pressure and Other Thermophysi- 
cal Properties 

The present determinations of TIT^ require measure- 
ments of the pressure and values for the thermal diffu- 
sivity of the gases; however, these quantities need not be 
known with nearly the same accuracy as the primary 
quantities. Expressions are presented for the virial coef- 
ficients, thermal conductivities, and viscosities of argon 
and xenon. 

7.1 Pressure Measurements 

The pressure was measured with a fused-quartz, 
bourdon-tube, differential pressure gauge (Ruska In- 
strument Corporation, Model No. 6000-801-1). The 
gauge had a full-scale range of 1 MPa and had a resolu- 
tion of 1 Pa. The reference side of the gauge was contin- 
ually evacuated with a mechanical pump to a pressure 
lower than 1 Pa. At the conclusion of the measurements 
the gauge was calibrated against an air-lubricated pres- 
sure balance designated "PG29" by the NIST Pressure 
Group (Ruska Instruments, Model 2460). The pressure 
indications from the bourdon-tube gauge were up to a 
factor of 4.5 X 10"* larger than the standard. The cali- 
bration data were represented by the equation 

/7/Pa = 7?b(l - 4.556 X lO""*), 

with a standard deviation of 4 Pa. (Here, R\, is the numer- 
ical value of the pressure indicated by the quartz 
bourdon gauge.) Between checks, the largest change of 
the zero-pressure indication of the gauge was 35 Pa; the 
average change between checks was 15 Pa. 

The gas in the resonator was separated from the man- 
ifold by a differential pressure transducer (DPT) con- 
structed of stainless steel and Inconel (MKS Instruments 
model 315BD-00100 sensor head with model 270B 
electronic display unit). The DPT had a full-scale range 
of ± 1 3 kPa and a resolution of better than 1 Pa. At zero 
differential pressure, the DPT's output voltage Vo varied 
with line pressure as 

yo/V = 2.744 X 10"X, 

where R^o is the numerical value of the pressure indicated 
by the quartz bourdon gauge. This equation was applica- 
ble up to 600 kPa and fit the results with a standard 
deviation of 1.2 mV (corresponding to 1.6 Pa). Between 
checks, the largest drift in Vo was equivalent to 12 Pa; 
however the average drift was 5 Pa. 

The DPT was calibrated over its full-scale range of 
±13 kPa at various line pressures between kPa and 



600 kPa before the acoustic measurements were made. 
The calibration depended upon both the sign of the 
differential pressure and the line pressure. The results 
for a positive differential pressure A+p were represented 
by the equations 

A+/7/Pa = (yc/V)[1.00213+ 2.299 X 10 '''7?b' 

- 2.469 X \Q-\Vjy)] 

with a standard deviation of 7 Pa. For negative differen- 
tial pressures A-/?, we found 

A-/7/Pa = (K/V)[0.99848 + 2.187 X \Q-'R^, 

+ 6.939 X 10-"(K/V)'] 

with a standard deviation of 7 Pa. In these equations, 
Vc=V — Vo is the voltage output of the capacitance 
differential pressure transducer corrected for the voltage 
output at zero pressure. 

We estimate that the standard uncertainty of a single 
measurement of the pressure of the argon was approxi- 
mately 19 Pa. For xenon, we checked the zero-pressure 
reading of the bourdon tube gauge after every pressure 
measurement. In this case, the standard uncertainty of a 
pressure reading was approximately 10 Pa. 

In this work, {T— T^o) was deduced by fitting polyno- 
mial functions of p to values of u^(p,T) and extrapolat- 
ing to u^(0,T) along isotherms. The pressure-dependent 
contribution to the uncertainty of (T—Tgo) appears on 
Row 13 of Table 2. For each isotherm, the contribution 
was estimated by multiplying the uncertainty of the zero 
of the pressure gauge by the linear term in u^(p ,T), i.e., 
Ai in Eq. (26). The pressure-dependent uncertainties are 
comparatively small. 

7.2 Thermal Conductivity and Viscosity of Argon 
and Xenon 

The thermal conductivity k of the gas affects the 
determination of T through the thermal-boundary-layer 
term in the corrections to the acoustic frequencies A/a. 
Since the uncertainties in the values of k at both T and 
Tw are correlated, the estimated 0.3 % relative standard 
uncertainty of k contributes only 0.2 X 10 "'' to the stan- 
dard uncertainty in T/T^, as indicated on Row 12 of 
Table 2. The viscosity 17 of the gas contributes to the 
half-widths of the acoustic resonances and to the very 
small terms in A/a resulting from crevices. Thus the 
uncertainty of the viscosity may be neglected when 
computing the standard uncertainty of T/T^. 
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We represented the thermal conductivity and viscos- 
ity of argon with the expression 

/<(r,p)/(mW-m-'-K-') = 6.4622 X 10""(r/K) 

- 7.7489 X 10-'(r/K)' + 5.4288 X lO'^iT/Ky 

+ 21.6X 10-'[p/(kg-m-')] (21) 

and 

'»7(r,p)/(|juPa-s) = 8.2822 X 10-"(r/K) - 1.0088 

X 10~'(r/K)' + 7.1206 X \0-'\T/Ky+ 11.10 

X 10-'[p/(kg-m-')] (22) 

respectively. For Eqs. (21) and (22), the zero-density 
values of k{T) and rjiT) were calculated from the HFD- 
B2 potential by Aziz and Slaman [24]. (The HFD-B2 
potential exploits theoretical results and simultaneously 
represents measurements of transport properties, den- 
sity virial coefficients, spectral data for argon dimers, 
and molecular beam scattering.) As discussed by Aziz 
and Slaman and in Refs. [4] and [7], the transport prop- 
erties derived from this potential are in agreement with 
recent measurements within experimental uncertainties, 
which range from 0.2 % to 0.3 %. The density coeffi- 
cients of the thermal conductivity and viscosity are the 
values tabulated by Maitland et al. [25]. 

The corresponding expressions for xenon are 

/<(r,p)/(mW-m-'-K-')= 1.8986 X 10""(r/K) 
- 1.9528 X 10-''(r/K)" + 6.0 X 10-'[p/(kg-m-')] (23) 
and 

'»7(r,p)/((xPa-s) = 0.079952(r/K) - 8.256 

X 10-''(r/K)" + 6.17 X 10-'[p/(kg-m-')] (24) 

respectively. For Eqs. (23) and (24), the zero-density 
values of k{T) and rjiT) were calculated from the HFD- 
B2 potential of Dham et al. [26] by Aziz [27] and the 
density coefficients were obtained from Ref. [25]. 

7.3 Density Virial Coefficients 

In order to correct the acoustic frequencies for the 
thermal boundary layer, we required the density p{T,p) 
and the constant-pressure molar heat capacity Cp{T,p) 
of the gases as functions of the temperature and the 
pressure. To obtain p{T,p), we inverted the virial equa- 



tion of state. We obtained Cp{T,p) from thermodynamic 
relations that use the virial coefficients and their tem- 
perature derivatives. 

For argon, the second virial coefficient B{T) and its 
temperature derivatives were calculated from the modi- 
fied HFD-B2 potential function of Ewing and Trusler, 
correct to the second quantum correction [28]. The 
parameters for this potential were determined solely 
from their recent precise measurements of the second 
acoustic virial coefficient of argon over a wide range of 
temperatures. Because of the underlying data, we expect 
that this potential function will return estimates of the 
virial coefficients that are more accurate than those 
from the FrFD-B2 potential of Aziz and Slaman [24]. 
The results were approximated by the polynomial 

5/(cm'-mor') = 34.236 - 1.166 X 10V(r/K) 

- 9.523 X 10'/(r/K)" 

with a negligible contribution to the uncertainties. Val- 
ues of third virial coefficient C were taken from [29] . 
The virial coefficients of xenon and its temperature 
derivatives were determined from the experimental data 
of Michels et al. [30]. The expression resulting from a 
polynomial fit to the data was 

5/(cm'-mor') = 58.509 - 0.3681 X 10V(r/K) 

- 0.5832 X 107(r/K)'. 

Values of C were also taken from Ref. [30]. 

8. Analysis of the Acoustic Data 

In this section we describe the analysis of the acoustic 
data. The resonance frequencies of the radial modes 
were corrected and then combined with values of ao{T), 
the radius at zero-pressure, that had been obtained from 
the microwave measurements. The resulting values of 
u(p,T)/ao{T) are listed in the Appendix in Table A2. 
The values of [u(p,T)/ao{T)f' were fitted to a polyno- 
mial (acoustic virial) expansion in the pressure to obtain 
the zero-pressure limits [uo{T)/ao{T)f'. The fitting was 
done in two, quite different, ways. First, the data for each 
isotherm were fitted by independent polynomials, and 
second, all of the data were fitted simultaneously to a 
u^(p,T) surface. The results of both fitting procedures 
were used to compute the ratios [uo{T)/ao{T)f/[uo{T^)/ 
«o(7"w)]^. To obtain T/T^, these ratios were combined 
with the radius-ratios determined from the microwave 
measurements and represented by Eq. (18). The temper- 
ature ratios were combined with the temperature mea- 
surements on ITS-90 to obtain (T—Tgo)- 
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The remainder of this section is subdivided as follows: 
Sec. 8. 1 describes the reduction of the raw data; Sec. 8.2 
describes the fits of the individual argon isotherms; Sec. 
8.3 describes the fitting of the argon u^{p,T) surface; 
Sec. 8.4 discusses variations of the fitting procedures, 
and Sec. 8.5 discusses the somewhat different analysis 
of the xenon data. 



8.1 Reduction of the Acoustic Data to Speeds of 
Sound on Isotherms 

The raw data for both argon and xenon were reduced 
to values of (u/a) on each isotherm in the following 
manner. The resonance frequencies were corrected for 
the thermal boundary layer, the effects of curvature of 
the surface, the temperature jump (assuming the ther- 
mal accommodation coefficient h=l), the coupling of 
gas and shell motions, and the effect of bulk dissipation. 
The resonance frequencies of the (0,2) mode were mul- 
tiplied by the factor 1 + 0.7 X 10"*^ to correct for the 
shape perturbation resulting from the unequal diameters 
of the hemispheres. All of these corrections are identical 
to those applied in previous work with this resonator and 
are given in Refs. [3] and [4]. As in Ref. [4], corrections 
were not applied for the small effects of the crevices 
surrounding the transducers because the geometry was 
not known well enough. The corrections for the static 
and dynamic compliances of the shell were taken from 
Ref. [2] and the dynamic compliance correction in- 
cluded the effect of radiation from the external surface. 
The elastic constants required for these corrections were 
taken from Ref. [4]; their temperature dependencies 
were taken from the measurements of Ledbetter et al. 
[31]. Finally, the frequencies were corrected to the exact 
isotherm temperatures T, using 

f{np)=f{T,p){TJTy" 

with negligible additional uncertainty. The corrected 
frequencies were divided by the appropriate eigenvalue 
and by UoiT). 

For tabulation only, the mean values of (w/ao) were 
computed for the first five radial modes. These means 
are listed in Table A2 of the Appendix along with pres- 
sures, temperatures, and the relative standard deviations 
of the mean. 

For fitting, the corrected frequencies were divided by 
the appropriate eigenvalue and by Co- At each tempera- 
ture and pressure there were two measurements made 
for each mode and the two values (m/ao) deduced from 
them were averaged. However, the values of (u/oo) for 
the five modes were not averaged so that the residual 
mode-dependence of the data could be studied. Follow- 
ing Refs. [4] and [7], each value of (u/ao) was weighted 



inversely by the square of its estimated standard devia- 
tion which was taken to be (o-j^ + 02^)"^. The first term 
(Ti is the estimate of the standard deviation of the fre- 
quency measurements from Ref. [7]: 

0-, = 1.4 X 10-'m'[1 + (10'Pa/;7)'(6 kHz//o„)']. (25) 

The second term, 02 = 3.7 X 10~'m^, accounts for the 
uncertainties of the mode-dependence of the model for 
the entire set of measurements. This term was chosen 
such that x^ = 1 for a "good" fit, as judged by deviation 
plots. The fitting procedures were adapted from Beving- 
ton [32]. 

8.2 Fitting Independent Acoustic Isotherms: 
Argon 

The corrected values of (u/aof were fitted by the 
expansion in powers of the pressure (in which we have 
omitted the subscript from the radius a for clarity) 

(u/af - (A,/a^)p^ = (Ao/a') + (AJa^)p + (Aj/a');?' 



+ {A-,/a^)p 



(26) 



from which we obtain Ao/a^ and, ultimately, T. 

In Eq. (26), the terms with positive integer powers of 
p are related to the virial coefficients that appear in the 
equation of state of a dilute gas. Thus, the form of these 
terms has a rigorous derivation, provided that the acous- 
tic frequency is sufficiently low, as is the case here, and 
provided the acoustic resonance frequencies are not ac- 
cidentally in near coincidence with a mechanical reso- 
nance of the steel shell of the resonator. (A possible 
coincidence of the (0,6) acoustic resonance with a me- 
chanical resonance is considered in Sec. 8.4.) In Eq. 
(26), the term iA-^ila^)lp relaxes the assumption that the 
thermal accommodation coefficient was exactly one. 
This assumption was part of the reduction of the reso- 
nance frequency data to speed-of-sound data on 
isotherms described in Sec. 8.1. 

Although the pressure on each isotherm was varied 
by a factor of 20, the values of (u/af on each isotherm 
varied by less than 0.3 %. At all the pressures used, the 
dominant term in Eq. (26), (Ao/a^), is far larger than the 
other terms. The smallest term, (A^p^la^), is always less 
than 3.1 X 10" ''of the dominant term (Ao/a^). Instead of 
trying to fit such a small term, we fixed the values of A3 
on the left hand side of Eq. (26) to results from Boyes' 
speed-of-sound measurements [33]. Boyes' measure- 
ments extended up to 10 MPa where (A^p^lA^^ becomes 
as large as 0.02. Boyes' data spanned the range 250 K 
to 350 K; thus, extrapolation to lower temperatures was 
required to determine A3 on our isotherms at 2 1 7 K and 
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234 K. To do so, we fitted a polynomial function of Tto 
Boyes' data. The resulting values of (A^/a^) are listed in 
the right-hand column of Table 7. We assumed that the 
uncertainty of the extrapolated A3 was equal to A3 itself. 
If this is true, the fractional effect of A3 on the values of 
Ao and on T deduced from Eq. (26) is only a few times 

For each argon isotherm, the parameters for Eq. (26) 
and the relative standard deviation of the fit are listed in 
Table 7. The mean of the relative standard deviations of 
the fits is 1.18 X 10 "^, which corresponds to approxi- 
mately 0.3 mK. The relative uncertainties of the 
parameters Ao/ao^ are carried to Rows 4 and 5 of Table 
2, where the uncertainties of (T—Tgo) are computed. 

Figure 9 displays the deviations of the data from Eq. 
(26) as a function of pressure. On the isotherms at 
273 K, 293 K, and 303 K, the deviations are clearly 
mode-dependent; the data for the (0,6) mode are farthest 
from the mean of the others, especially at the higher 
pressures. This will be discussed in Sec. 8.4. 

8.3 Fitting M^(/;,r) Surface 

The isotherm analysis is one extreme, insofar as it 
uses the fewest possible assumptions concerning the 
fitting parameters. Thus, the results do not have uncer- 
tainties resulting from choosing only one among several 
competing models for the data. Here, we consider an 
alternative that tends towards the opposite extreme: we 
impose as many physically-based constraints on the 
parameters as possible. Remarkably, the imposition of 
these constraints changes (T—Tgo) on the average by 
0.6 mK, and 0.8 mK in the worst case. Furthermore, 
these changes are within the combined uncertainties of 
the two analyses (see Table 1). We expect that imposing 
more constraints than in Sec. 8.2 but fewer constraints 
than imposed here will lead to results between those 
resulting from these extremes. 

The most significant constraint when fitting u^(p,T) 
was applied to At{T). We noted that At{T) = y°l3^{T)/M 
where I3^{T) is the second acoustic virial coefficient. 
To = Cp°/Cy° = 5/3 is the ideal-gas heat-capacity ratio for 



argon, and M is the molar mass of the argon. One may 
calculate fS^T) from the argon-argon interatomic pair 
potential. The semi-empirical HFD-B2 potential func- 
tion of Ref [24] has been refined to simultaneously 



< 



I 



n 



J O o 



. « 



tl 



I ' ' 

° 9 « 





• 


(0,2) 





(0,3) 


T 


(0,4) 


V 


(0,5) 


■ 


(0,6) 



r=217K 



^ V ▼ 

' 8 • o 
9 5 



r=234K 



5 o 



S ■ - ■ ® * Q 

1 5 9 § i » ! 



r=253K 



♦ n « 0» 



e 



r=273K 



_ ^ ■ 



- V • 



1999 



r=293K 



8 • 8 



§ ? § 



9 
r=303K 



100 200 300 400 500 600 

p/kPa 

Fig. 9. Speed of sound in argon: isotherm analysis. Deviations of the 
data for the speed of sound in argon from Eq. (26) with the parameters 
from Table 7 (A«^ = unmeasured - «\tted). 



Table 7. Parameters from analysis of the acoustic isotherms of argon 



T 


Ao/ao 


lO'^AiW 


IO'AjW 


lO^M-iW 


10"^ (t{u^)Iuo^ 


10"*A3W 


(K) 


(s-') 


(s^^-Pa"') 


(s-n-Pa-^) 


(s-n-Pa) 




(s^^-Pa"') 


302.9166 


13 282 850 ± 19 


66 121 ± 95 


6.01 ± 0.12 


76 ±94 


1.28 


60 


293.1300 


12 857 689 ± 16 


54 820 ± 83 


6.25 ±0.11 


74 ±76 


1.03 


87 


273.1600 


11 989 200 ±20 


28 513 ± 109 


6.73 ± 0.16 


18 ±94 


1.32 


147 


253.1500 


11 117 708 ± 15 


-2 868 ± 83 


6.87 ± 0.16 


84 ±65 


0.88 


216 


234.3156 


10 296 361 ± 19 


-38 619 ±98 


6.64 ± 0.14 


58 ±85 


1.19 


293 


217.0950 


9 544 464 ± 22 


-78 398 ± 150 


6.27 ± 0.29 


59 ± 89 


1.21 


374 
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represent information from theory, molecular beam 
scattering, spectroscopy of argon dimers, measure- 
ments of density virial coefficients, and measurements 
of transport properties. However, it does not make use 
of speed-of-sound data. Ewing et al. measured fS^T) for 
argon with an acoustic resonator [34]. Remarkably, their 
values of fS^T) differed from those computed from the 
HFD-B2 potential function by a very small, linear func- 
tion of the temperature in the range from 100 K to 
304 K. We used Ewing et al.'s modified HFD-B2 poten- 
tial function to compute the values A,{T) that served as 
a fixed base line when fitting the present data (these 
values appear in Table 8). When fitting our data, we 
included the term 



AA,(T)/ao^^bo + b,(T-TJ. 



(27) 



Thus, we constrained the values of Ai{T) on six 
isotherms to be fitted by two parameters, bo and bi, 
instead of six parameters. 

Table 8. Parameters from the u^(p,T) surface analysis for argon. The 
thermal accommodation coefficient h was constrained to be exactly 1 . 
A\{T) differs from Ai(r) predicted by the HFD-B2 interatomic poten- 
tial by a fitted linear function [Eq. (27)] of T. Also, A2{T) is fitted by 
a quadratic function [Eq. (28)] of l/T. 



T 


Ao/ao' 




lO'AiW 


10 


'"'As/ao' 


(K) 


(s-^) 




(s"^-Pa"') 


(s 


-^■Pa-^) 


302.9166 


13 282 869 ± 5 




64 879 




60 


293.1300 


12 857 707 ± 4 




53 520 




87 


273.1600 


11 989 191 ± 3 




27 191 




147 


253.1500 


11 117 724 ±4 




-4 437 




216 


234.3156 


10 296 389 ± 4 




-40 442 




293 


217.0950 


9 544 478 ± 5 




-80 295 




374 




10X = (13.67 ± 


0.23) 


s-'-Pa-' 







lO^fo, = (-7.31 ± 0.77) s" -Pa" -K"' 
10% = (-15.31 ± 0.90) s"^-Pa"^ 
lO'^c, = (10.91 ± 0.46) s"^-K-Pa"^ 
10'c2 = (-1.339 ± 0.061) s"^-K^-Pa" 



The second constraint when fitting u^(p,T) was ap- 
plied to A2iT). We used the expression 



A2{T)/ao^ = Co + Ct/T+ cJT^. 



(28) 



thereby constraining the values of A2iT) on six 
isotherms to be fitted by three parameters: Co, c-i, and C2. 
As mentioned in Sec. 1, the rationale for this constraint 
begins with the observation that the present isotherms 
are well above the critical temperature T^ of argon 
(1.4 < T/T^ < 2.0) where the viral coefficients of argon 
are only weakly temperature dependent. Furthermore, 
the present data are well below the critical density 
(p/p^ < 0.02), where the virial expansion converges 



rapidly. The ratio of terms {A2p^la^)l{AJa^) attains a 
maximum value of only 1.5 X lO"'*; thus, only a moder- 
ately precise representation of A2{T) is required. 

The final constraint for fitting the u^(p,T) surface 
was to assume that the thermal accommodation coeffi- 
cient h is exactly 1, which is its temperature-independent 
maximum value. (If we had simply assumed that h was 
independent of temperature, its optimum value would 
have been 0.90 and the resulting values of T/T^ would 
have been negligibly different from the present surface 
fit.) 

The parameters resulting from fitting the u^(p ,T) sur- 
face to the argon data are listed in Table 8. The fractional 
uncertainties of the parameters Asja^ are carried to 
Rows 6 and 7 of Table 2, where the uncertainties of 
{J— Tgo) are computed. The deviations of the data from 
the fit are shown in Fig. 10. The deviations are compara- 
ble to those obtained from fitting the data on each 
isotherm separately (See Fig. 9.). The relative standard 
deviation of the fit was 1.12 X 10"^, which is close to 
1.18 X 10"^, the mean of the relative standard devia- 
tions from the isotherm fits. Thus, by statistical criteria, 
the 1 1-parameter surface fit to the u^(p ,T) data is just as 
valid as the isotherm fits that require 24 parameters. 

8.4 Alternative Fits of Argon u^(p,T) Data 

As evident in Fig. 9 for the isotherms at 273 K, 
293 K, and 303 K, the results from the (0,6) acoustic 
mode differ from the average results from the other 
modes. A similar phenomenon was observed during the 
previous redetermination of T^ with this resonator [7]. 
Then, the (0,6) mode was perturbed more strongly at 
10.0 kHz at Tg than at 9.5 kHz at r„ and the perturbation 
was attributed to the near coincidence of the (0,6) fre- 
quency with a non-radial shell resonance that had been 
predicted to occur at 10.15 kHz. Perhaps the intervening 
modifications of the apparatus, such as attaching the 
copper shield to the resonator, reduced the frequency of 
the non-radial shell mode. If this had happened, the near 
coincidence of the shell mode and the (0,6) mode might 
occur at the lower temperatures where the frequency of 
the (0,6) mode was also lower. 

In Ref. [7] an additional parameter, A/, applicable to 
the (0,6) mode only, was included in the fit to account 
for the coincidence of that mode with the shell reso- 
nance. In this work, by contrast, we elected to repeat 
both the isotherm and the surface analyses of the acous- 
tic data excluding the (0,6) data entirely. Table 9 com- 
pares results from these analyses with the earlier ones. 
Upon excluding the (0,6) mode, the values of (T—Tgo) 
are reduced only about 0. 1 mK and the standard devia- 
tions of the fits are reduced by approximately 25 %. 
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Fig. 10. Speed-of- sound in argon: surface analysis. Deviations (A«^ = «^„eiisured ~ "^aicd) of 
the data from Eqs. (26), (27), and (28) with the parameters from Table 8. The results for the 
different radial modes are denoted (0,2) etc. 



Table 9. The difference T-T. 



90 







Surface fits 




Isotherm fits 




Sec. 
8.3 


(0,6) 
omitted 


quadratic 

AAi(r) 


Sec. 8.2 


(0,6) 
omitted 


10''(7(hV«o' 


1.12 


0.84 


1.11 


1.18 


0.85 


Tso/K 






T-T^ 







302.9166 
293.1300 
253.1500 
234.3156 
217.0950 



4.61 

3.23 

-2.29 

-2.92 

-3.55 



4.48 

3.10 

-2.43 

-3.06 

-3.65 



4.25 

3.06 

-2.34 

-3.12 

-3.92 



3.95 

2.60 

-2.84 

-3.73 

-4.03 



3.82 

2.56 

-2.89 

-3.72 
-4.23 



For evaluating this result and for future work, it is 
useful to compare the effect of the compliance of the 
gas-constant resonator's shell to the coefficient Aj. This 
was done in Eq. (3.1) of Ref. [4] for the radially-sym- 
metric ("breathing") mode of the shell. A radial acous- 
tic mode at frequency /a is shifted by A/shdi, where 



A/sh 



-roA's.ip/3 



/a 1 - (/a/13.58kHz)^ 



(29) 



and where the compliance for internal pressure is 
A-M = (3.03 ± 0.03) X 10"" Pa"' at r„. Far from 
13.58 kHz, A/shell is a linear function of pressure that 
contributes 17 X 10"*^ to AmVm^ near 500 kPa. If A/sheii 
were not accounted for, its effects would be obvious on 
deviation plots, such as Fig. (9). Upon comparing A/sheii 
to Ai at, for example, r„ we find ('yoA's,i/3)Mi ~ 0.007. 



Thus, far from the resonance, the effect of radial com- 
pliance is twice the uncertainty of Ai as deduced from 
the isotherm fits and close to a resonance the effect will 
be much larger. Similar considerations apply to the 
bending modes of the shell that fall below 13.58 kHz. In 
cases where a shell resonance is suspected, the reso- 
nance could be unambiguously revealed by the acquisi- 
tion of data on closely spaced isotherms. In the vicinity 
of the resonance, A\(T) will have an anomalous temper- 
ature dependence for the acoustic mode closest to the 
shell mode. When /a is very close to a shell resonance, 
expressions similar to Eq. (29) must be superceded by 
more complex ones that account for damping of the 
shell's motion and for the "avoided crossing" of the gas' 
resonance and the shell's resonance. 

We also considered a surface fit such that AAi(r) had 
the quadratic representation 



36 



Volume 104, Number 1, January-February 1999 

Journal of Research of the National Institute of Standards and Technology 



^A,{T)lal ^bo + b,iT-TJ + b^iT-T^f 



(30) 



The results fell approximately half way between the 
surface fit and the isotherm fits with a negligible reduc- 
tion in the standard deviation of the fit (see Table 9). 

The fits that we considered here and the analysis of all 
the previous results from this resonator weighted the 
acoustic data with the pressure dependence given by 
((Ti^ + a-2^)~\ where (Tj and a-2 are given by Eq. (25) and 
the text associated with that equation. However, the sig- 
nal-to-noise ratio of the present data does not decrease 
with pressure according to Eq. (25) and previous experi- 
ence with this apparatus [see Fig. (10)]. Presumably, the 
improvements mentioned in Sec. 4 account for the 
change. In future work, we will investigate the effects of 
weighting the acoustic data independently of the pres- 
sure and we will implement the improvements enumer- 
ated in Sec. 1.1. We hope that these changes in the 
apparatus, the procedures, and the weighting of the data 
will determine the parameter {AJa^) in Eq. (26) more 
precisely. Perhaps these changes will reduce the differ- 
ences between the surface fits and the isotherm fits in 
Table 9. If so, the overall uncertainties in (T—Tsio) might 
be reduced. 



8.5 Acoustic Data for Xenon 

8.5.1 The Virtual Leak 

The most important problem encountered when ana- 
lyzing the acoustic data for xenon was the progressive 
change in the speed of sound in the gas while the gas 
was in the resonator. Because there were no known 
leaks in the gas handling system, we attributed the con- 
tamination to a virtual leak, that is, a volume sealed from 
the laboratory and connected to the resonator by a path 
of low pumping speed. If such a volume were exposed 
to a contaminating gas at high pressure, it would fill 
rapidly via Poiseuille flow. Subsequently, such a volume 
would be a lingering source of contamination because it 
would evacuate comparatively slowly via molecular 
flow. A precursor to this problem was detected during 
the redetermination of R. We quote from Ref. [4]: 
"When the resonator was filled with helium, we mea- 
sured a slow decrease in the resonance frequencies. In a 
typical case the fractional decrease was 9.3 X 10~V 
(100 h) with 438 kPa of He-M in the resonator. . . We 
speculate that slow desorption of impurities is responsi- 
ble for these effects. Possible sources of water etc. are 
the "Viton" 0-rings which seal the microphone ports 
and the fill port to the resonator." To this list of possible 
sources, we add the wax used to seal the two hemi- 
spheres together. 



During the present measurements, when the res- 
onator was filled with xenon at Tg, the resonance fre- 
quencies slowly increased. In a typical case, the frac- 
tional increase was 2 X 10 V(60 h) at T^ with 150 kPa 
of Xe in the resonator. Presumably, the xenon was pro- 
gressively contaminated with an impurity with a higher 
speed of sound. Again, we suspected that the impurity 
had been dissolved in the 0-rings. Perhaps, for example, 
argon had diffused into the 0-rings comparatively 
rapidly while the resonator was filled with argon at high 
pressure. Later, when the argon in the resonator was 
replaced by xenon, small amounts of argon slowly des- 
orbed from the 0-rings into the resonator, driven only 
by the concentration gradient remaining within the O- 
rings. This speculation is consistent with measurements 
conducted with argon immediately following the xenon 
measurements. Then, a progressive fractional decrease 
in the resonance frequencies of 4 X 10 'V(40 h) was 
measured at T^ with 150 kPa of argon in the resonator. 
The change in sign is consistent with the argon sample 
being progressively contaminated by desorbing xenon 
and the larger rate is reasonable considering that 
{l/u^){du^/dx) is 3.3 times larger for dilute xenon in 
argon than for dilute argon in xenon. In contrast, pro- 
gressive contamination was not detected during similar 
measurements at r„ with either argon or xenon in the 
resonator. We do not know whether the capacity of the 
0-rings is much lower at T^ and/or the rate of desorption 
is so slow (or so rapid) that we did not observe it. 

We emphasize that when the resonator was filled with 
argon prior to the xenon measurements, we never ob- 
served a secular change in the frequencies at any of the 
temperatures under study. 

8.5.2 Correction of the Xenon Data for 
Contamination 

In order to cope with the progressive contamination of 
the xenon at T^, we decided to measure the resonance 
frequencies as a function of time t at each pressure. The 
values of /a(f) were extrapolated backward to the time 
when the resonator was first filled and when, pre- 
sumably, the contamination began. To minimize the un- 
certainty resulting from the extrapolation, f^{t) was 
measured at T^ for three separate fillings of the res- 
onator. This greatly reduced the residence time of each 
sample in the resonator and hence its progressive con- 
tamination. Furthermore, the resonator was evacuated 
and "baked" between successive fillings. 

At each temperature and pressure, the values of the 
derived speeds of sound from each mode were averaged. 
The averages were corrected for progressive contamina- 
tion assuming that, for each sample, the contamination 
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effect was inversely proportional to the pressure and 
directly proportional to the sample's residence time at 
each pressure. For each sample, successive corrections 
for the different pressures were assumed to be additive 
and the time required to reduce the pressure was ne- 
glected in comparison with the residence time at any 
particular pressure. Thus, the correction Am for the (-th 
pressure is 



Am,- = Y,ktjlpj 



(31) 



j=i 



in which tj is the time spent at pressure pj . We did not 
assume that the outgassing-rate-constant k was the same 
for each sample. Instead, we assumed that it was inde- 
pendent of pressure and time for a given sample and we 
determined it by monitoring the drift of the resonance 
frequencies at constant temperature and at a low pres- 
sure for at least 1 2 h. The values of k for each of the 
three samples are listed in the right column of Table 10. 
Table 10 also contains the averaged values of u{p,T^,t) 
and the values extrapolated to t = 0. 

The extrapolation procedure was successful. The de- 
viations of the xenon data at T^ from fitting functions 
do not show discontinuities where the samples were 
changed (see below and Fig. 11). The maximum frac- 



tional correction to 



for the virtual leak was 



7.6 X 10"'' with an uncertainty of approximately 
0.24 X 10"^. If the correction had been neglected, the 
derived value of T„/T^ would have been erroneous be- 



cause of an incorrect determination of the curvature of 
the isotherm. 

8.5.3 Fitting to Xenon Acoustic Data 

The xenon data for the isotherms at r„ and T^ were 
corrected for all of the phenomena that were discussed 
in Sec. 8.1 in connection with fitting the argon data for 
the six isotherms. Then the xenon data were averaged 
and corrected for the virtual leak. The resulting values 
of u/oq averaged over five modes are listed in Table A3 
in the Appendix. 
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Fig. 11. Speed-of-sound in xenon. Deviations (A"^ — ^^mcasmd ~ 
"^fitted) of the data from Eq. (26) with A-i = and with the other 
parameters from Table 1 1 . 



Table 10. Virtual leak correction for xenon at 302.91 K 



Run 


P 


t 


(«/«>, 


10'^8« 


(«/a>,_,o 


lOV 


10*^*: 


no. 


(kPa) 


(h) 


(s-') 






(s-') 




(kPa-h"') 


1 


33.095 
33.065 


5. 
19.5 


2009.2915 
2009.2952 


1.46± 


0.64 


2009.2912 


0.86 


3.32 


2 


302.569 


16. 


1998.9358 






1998.9355 


0.02 






269.851 


19. 


2000.1993 






2000.1990 


0.03 






239.836 


22. 


2001.3573 






2001.3569 


0.04 






199.948 


25. 


2002.8936 






2002.8931 


0.05 






150.100 


28.5 


2004.8097 


0.95 ± 


0.22 


2004.8092 


0.06 


2.39 




150.038 


88. 


2004.8141 












3 


133.515 


13.5 


2005.4467 






2005.4456 


0.13 






105.227 


17. 


2006.5318 






2006.5302 


0.17 






79.949 


20.5 


2007.5001 






2007.4980 


0.23 






54.958 


23.8 


2008.4570 


1.36± 


0.30 


2008.4542 


0.29 


5.75 




54.941 


36.5 


2008.4600 














44.964 


40.8 


2008.8429 






2008.8364 


0.62 
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Equation (26) was independently fitted to the cor- 
rected data for each isotherm, as discussed in Sec. 8.2, 
except that the characteristic pressure in the weighting 
function [see Eq. (25)] was replaced by 5.5 X lO"* Pa 
and the value of A3 was fixed at zero. (The maximum 
value of {A^/Ao)p^ is only 9 X 10~^ in the present work. 
This was determined by estimating A3 for xenon from a 
dimensionless plot of (A^/AoXTpJT^y as a function of 
T/T^ that we made with acoustically-determined values 
of A3 for argon [33] and xenon [35].) 

Because there were only two xenon isotherms, we 
considered only one constraint; namely, fixing the ther- 
mal accommodation coefficient at unity or, equivalently, 
requiring A. 1 to be identically 0. Upon allowing A-i ¥= 
0, the standard deviation of the fits increased and the 
best values of A_i did not change from zero within their 
uncertainties. Thus, both the xenon data and the argon 
data are consistent with A-i identically equal to 0. 

The parameters and the standard deviations resulting 
from the fits are listed in Table 1 1 . Figure 1 1 shows the 
deviations of the xenon data from Eq. (26) for the con- 
strained fits at Tw and T^. For the constrained fit, the 
uncertainties in AJa^ contribute 1.3 X 10"'' to the un- 
certainty in TJT^; for the unconstrained fit, the contri- 
bution is 4.6 X 10"^. When the non-acoustic sources of 
uncertainty (Table 2) are added in quadrature, the con- 
strained xenon fit leads to the result (Fg— Tgo) = 4.38 ± 
0.48 mK and the unconstrained fit leads to (T^—Tsio) = 
4.38 ± 1.42 mK. 

Most of the deviations in Fig. 1 1 are independent of 
pressure and they do not show discontinuities at the 
pressures where the resonator was refilled. The devia- 
tions do depend upon the mode index in a similar man- 
ner on both isotherms. This led us to consider small 
adjustments to the eigenvalues that could be rationalized 
by our imperfect knowledge of the geometry of the 
resonator. Such adjustments reduced the standard devia- 
tion of the xenon fits by a large factor. However, the 
same adjustments had essentially no effect on the argon 
surface fit. We concluded that the xenon deviations are 
a measure of our ignorance of the resonator's behavior. 
It may be useful to reconsider this matter in the future, 
when argon data become available on closely spaced 



isotherms. Such data may separate the effects of imper- 
fect knowledge of the eigenvalues from imperfect 
knowledge of the shell's dynamic compliance. 



9. Comparisons with Previous Acoustic 
Measurements 

9.1 Comparison at r„ 

In the upper panel of Fig. 1 2, the present results at T^ 
are compared with previous measurements made with 
the same gas in the same apparatus. The baseline is the 
fit to the isotherm used to redetermine R and reported in 
Table 1 1 of Ref. [4]. Obviously, the differences between 
the present results and the measurements conducted in 
1986 are very small; the average difference is on the 
order of, fractionally, 1 X 10"'' over the entire pressure 
range. It is possible that some of this difference could be 
explained by recalling that, following the 1 986 measure- 
ments, a transducer was relocated resulting in a possible 
decrease in the volume of the resonator by as much as 
0.88 X 10"'^, fractionally (see Sec. 9.6 of Ref. [4]). This 
corresponds to a 0.59 X 10"'' fractional increase in u^. If 
so, the discrepancy is within the experimental uncer- 
tainty. 

The upper panel of Fig. 12 also displays the unpub- 
lished results of 1989 [8]. They are close to both the 
present results and those of 1986. This figure provides 
convincing evidence that the volume of the gas-constant 
resonator has been very stable for many years. 

9.2 Comparison at T^ 

The lower panel of Fig. 12 compares the present re- 
sults for M ^ at Tg with the earlier ones obtained with the 
same apparatus. For this panel, the baseline is Eq. (26) 
with the present parameters from Table 8. The present 
measurements nearly agree with the unpublished results 
of 1989. The maximum fractional difference is only 
2 X 10"^ at the highest pressure. This difference de- 
creases approximately linearly with pressure, and the 
zero-pressure values of m ^ differ by less than 1 X 1 0"^ m ^, 



Table 11. Parameters from analysis of the acoustic isotherms of xenon 



T 


Ao/oo^ 


lO'^AiW 


IO'AjW 


lO^^A-iW 


10' 


a{u')luo 


(K) 


(s-^) 


(s-^-Pa"') 


(s-'-Pa-^) 


(s-^-Pa) 






273.1600 


3 648 625.6±3.3 


-204 000 ±39 


-8.91±0.10 







1.28 


302.9166 


4 042 322.7 ±3. 7 


-153 363 ±46 


-1.93±0.13 







1.40 


273.1600 


3 648 617.1±11.5 


-203 940±86 


-9.03±0.18 


31±39 




1.28 


302.9166 


4 042 319.8±14.1 


-153 342±111 


-1.98 ±0.24 


10±47 




1.51 
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Fig. 12. Comparison of present results with previous results obtained 
with the same apparatus (Am^ = M^„easurcd ~ wetted)- Top: u^aud for this 
plot came from the redetermination of R [4]. Bottom: u^aud is Eq. (26) 
with the parameters for Tg from Table 8. The "Ref. [7] 1988" data were 
used to remeasure Tg and are now suspect. 

which is equivalent to only 0.3 mK in (T—Tgo)- This is 
remarkably good agreement. 

In contrast, the data labeled "Ref [7] 1988" from the 
previous redetermination of T^ are seen to be in poor 
agreement with the present results. At the highest pres- 
sure the results agree to better than 2 X 10"'', fraction- 
ally. However, as the pressure decreases, the difference 
grows non-linearly, reaching approximately 9 X lO"'', 
fractionally, at 25 kPa. We now believe that this was a 
result of progressive contamination of the argon sample, 
similar to the contamination of the present xenon sam- 
ples and the earlier helium samples that we described in 
Sec. 8.4.1. 

One question remains. Common impurities increase 
the speed of sound in argon; however, the departure of 
the earlier data from the present data is consistent with 
a progressive decrease of the speed of sound. The ques- 
tion is: what impurity could have caused this decrease? 

We considered the possibility that a droplet of mer- 
cury remained in the resonator from the weighings asso- 
ciated with the redetermination of R. The saturated 
vapor pressure of mercury is 0.37 Pa at Tg and 
(l/M^)(dM^/dx) = — 4 for a mercury impurity in argon. 
Thus, mercury at its saturated vapor pressure would 
fractionally decrease the speed of sound in argon at 
100 kPa by 15 X 10 "^, much more than was observed. 



We note from Sec. 9. 1 of Ref. [7] that helium was in 
the resonator prior to the argon measurements in ques- 
tion and that the helium measurements were followed by 
baking out the resonator under vacuum for 4 d at 60 °C 
and two flushes with argon. This probably was the first 
time that the gas constant resonator had been maintained 
at such a high temperature since it was assembled. Per- 
haps some moderately heavy hydrocarbons diffused out 
of the wax that sealed the two hemispheres together and 
dissolved into the O-rings. The experience with xenon 
indicates that flushing did not remove the source of the 
contamination. Perhaps the same phenomena occurred 
during the "Ref [7] 1988" argon measurements. 



10. Checks for Systematic Errors 

In this section we consider tests for systematic errors. 
In the first test, the isotherm parameters determined in 
this work are compared with previously published val- 
ues. In the second, the measured half- widths are com- 
pared with those calculated from the theory. In the third 
test, the molar mass of xenon calculated from the 
present results is compared with the value determined 
from lUPAC tables of the isotopic composition of xenon. 
In all cases we find satisfactory agreement. 

10.1 Discussion of the Isotherm Parameters 

10.1.1 The Regression Coefficient Ai{T) 

For comparisons, we consider the more familiar sec- 
ond acoustic virial coefficient ^JJ") = RTA\{T)IAa{T). In 
Fig. 13, we compare the values of A;Sa = [^JDm^^nKi - 
;Sa(r)caicuiated] from our Independent isotherm analysis 
(Table 8) and from other measurements. The baseline of 
Fig. 1 3 is the values of ^JJ) that were calculated from 
the HFD-B2 intermolecular potential of Ewing and 
Trusler [34]. This potential was adapted by Ewing et al. 
[36] to fit their own measurements of ^JJ) that are 
labeled "UCL 1989" on the figure. The uncertainty of 
the baseline must be on the order of the experimental 
uncertainty associated with their results (rms uncer- 
tainty of 0.081 cm^-mol"'). 

In the present work, the largest contribution to the 
uncertainty of ^JJ) is the standard deviation of Ai(r) 
from the isotherm analysis. The standard deviations in 
Table 7 propagate into uncertainties of 0.010 cm^-mol~' 
to 0.022 cm^-mol~' for ^JJ). The present values of 
PJJ) are subject to two sources of uncertainty arising 
from systematic effects that were discussed in connec- 
tion with the redetermination of R [4]. First, the uncer- 
tainty of the shell's compliance ;^'s,i was estimated to be 
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Fig. 13. Deviation of experimental values of Pa{T) from the values 
calculated from the potential of Ref. [34]. Key: •, this work; X, Ref. 
[4]; +, Ref. [7]; V, Ref. [37]; A, Ref. [33]; O, Ref [35]; D, Ref. [35]; 
A, Ref [38]. 

6 % of its value, independent of the temperature. This 
propagates into an uncertainty of 0.085 % ofAi(T^) and 
it is equivalent to an uncertainty of 0.0046 cm^-mor ' in 
I3JJ"). Second, the uncertainty of the compliance of the 
transducers propagates into an uncertainty of 0.0025 
cm^-mol~^ in I3JJ"). We calculated the standard uncer- 
tainty of PJT) from the sum in quadrature of these two 
terms and the standard deviation of Ai(r) from the fit. 
We also included terms from the uncertainties of h and 
of the thermophysical properties required in the calcula- 
tions; these contribute less than 0.01 cm^-mol~' to the 
uncertainty of I3JJ"). 

Figure 13 shows that the present results differ from 
the potential of Ewing and Trusler [34] by approxi- 
mately 0.09 cm^-mol~' across the whole range. This is 
barely outside the combined uncertainties. The value 
determined at T^ differs from the value from [4] by only 
0.009 cm^-mor\ This level of agreement is truly re- 
markable. In contrast, this determination of ;Sa at T^ 
differs from that of Ref. [7] by 0.18 cm'-mor\ The 
results determined from the high pressure measure- 
ments of Ewing and Goodwin [37] are shown together 
with results from Goodwin and IVIoldover [38], and from 
Boyes [33]. The agreement is better than 0.2 cm^-mol~^ 
over the whole range. This level of agreement is extraor- 
dinary considering the pressure ranges under study dif- 
fer by more than a factor of 10. 

10.1.2 The Regression Coefficient AiiT) 

For comparisons, we converted the present values of 
A2(T) to the third acoustic virial coefficient yJ.T) = 
RTA2(T)/Ao(T). In Fig. 14, we compare the values of 
yJ.T) determined in this work with previous values. The 



value of yJ.T^) measured when R was redetermined [4] 
agrees with the present value to within the experimental 
uncertainty. The results of Ewing and Goodwin [37], 
and of Boyes [33], which extend to pressures far greater 
than were studied here, agree with the present results to 
within twice the combined experimental uncertainties. 
However, the analysis of the present data incorporated 
Boyes' values of A3; thus, the measurements are not 
completely independent. The value of yJ,T^ from [7] 
labeled "NBS 1988" is not consistent with the other 
values of yJ,T). This is another indication that the 
IVIoldover-Trusler redetermination of T„ was in error. 







1 




1 

r 


1 1 
1 

A 

9 


1 


1 


1 1 




1.30 

"0 

E 


1 
~ 41 


V 




• This worl< 
X NBS 1987 
+ NBS 1988 
VUCL1988 
AUCL1992 




T 1.20 
re 

Q- 

_E 

^^° 1.10 




A 

















1.00 


" 








A - 




1 


1 1 








250 



300 



350 



r/K 



Fig. 14. Comparison of values of yaiT). Key: •, this work; X, Ref. 
[4]; +, Ref. [7]; V, Ref. [37]; A, Ref [33]. 



10.2 Half- Widths of the Resonances 

At low pressures, the largest correction to the mea- 
sured resonance frequencies is determined by the acous- 
tic boundary layer and this same boundary layer is also 
the largest contributor to the half-widths of the acoustic 
resonances. Thus, the measurements of the half -widths 
provide an important method for confirming the appli- 
cability of the theory of the boundary layer to the 
present determination of T/T^. We define the excess 
half-width Ag as the amount by which the measured 
half-width exceeds the calculated half -width. Figure 15 
displays Ag for five radial modes in argon at r„ and T^. 
(The values of Ag are scaled by lOV/a, where /^ is the 
resonance frequency of the mode under study.) For each 
mode on Fig. 15 , an extrapolation top = yields values 
of Ag/f< lO"''. This is strong evidence that the 
boundary-layer correction is applicable to the present 
data at the lower pressures. The results on Fig. 15 are 
similar to those observed previously with either argon or 
helium in this resonator at T^ and at T^ and similar to the 
results obtained in this work on the isotherms below T^. 
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Fig. 15. Scaled excess half-widths of the resonances in argon as a 
function of the pressure (Ag = g„„sured - gcaicukted)- 

Figure 16 displays Ag// for the measurements with 
xenon at r„ and Tg. At low pressures, Ag//of the (0,7) 
mode increases sharply as a result of partial overlap 
with the neighboring non-radial (13,2) multiplet (see 
Fig. 8 of Ref [4]). The data for the (0,8) mode also 
indicate a problem upon extrapolation to /? = 0. For the 
five lower modes, Ag//< 2 X 10 "^ at the lower pres- 
sures. This bound is small enough to indicate that the 
boundary-layer correction for the xenon data is not so 
seriously in error as to vitiate the current redetermina- 
tion of Tg/r^ with xenon. However, the pressure depen- 
dence of Ag// is more complicated than that for argon 
and may be a worthy subject for further study, especially 
under circumstances where progressive contamination is 
not an issue. 



10.3 Molar Mass of Xenon 

The regression parameter Ao = u^ = y°RTIM was de- 
termined at both r„ and Tg for both argon and xenon. 
From these values it is possible to determine the ratio of 
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Fig. 16. Scaled excess half-widths of the resonances in xenon as a 
function of the pressure (Ag = g„„sured - gcaicukted). 

the molar mass of our working xenon sample M„xe to 
the molar mass of our working argon sample M„Ar- We 
have two estimates of M^xJM^at and the difference be- 
tween them is an indication of uncertainty, probably 
resulting from the contamination of the xenon. We find 
M„xe/M„Ar = 3.285 957 ±0.000 004 from the mea- 
surements at r„; we find M„xe/M„Ar = 3.285 949 ± 
0.000 005 from the measurements at T^. The fractional 
difference between these ratios is 2.4 X 10"^ which of 
course is the fractional difference between the values of 
TJT^ determined with xenon and with argon. 

The molar mass of our argon sample was determined 
previously [4] with a relative uncertainty of only 
0.8 X lO"'' by comparison with a sample of chemically 
purified, nearly mono-isotopic ''"Ar. Using the value of 
(23.968 684 ± 0.000 019) g/mol for M„Jy° we find 
M„xe//= (78.759 98 ±0.000 23) g/mol, where we 
have averaged the two ratios. Our sample of xenon con- 
tained krypton with a concentration of 18 X lO"'', by 
mole fraction. Accounting for this leads to a value for the 
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molar mass of Xe of (131.267 49 X 0.000 33) g/mol 
where we have taken y" = 5/3 exactly. Within com- 
bined uncertainties, this value agrees with the value of 
(131.29 ± 0.03) g/mol recommended by lUPAC. Ref- 
erence [39] states that "modified isotopic compositions 
may be found in commercial material ... because the 



material has ... inadvertently been subjected to isotopic 
separation" leading to systematic differences in the esti- 
mated molar masses of various different samples. On 
the scale of the lUPAC uncertainty, our problems asso- 
ciated with contamination of the xenon are minor in- 
deed. 



11. Appendix A. Tabulated Data 



Table Al. Frequencies of microwave resonances 



TM/m 


/d/MHZ 


/,/MHz 


(/"„>/MHz 


gm(expt) 
(MHz) 


gm(calc) 
(MHz) 


10'(Ag//) 


lO'ecjo 


10''ba/a 








T/K = 


= 210.7400 










11 
12 
13 


1473.6870 
3285.6832 
5004.5941 


1474.0150 
3286.0910 
5005.1579 


1473.7963 
3285.8191 
5004.7820 


0.1392 
0.1594 
0.1909 


0.1233 
0.1428 
0.1708 


10.8 
5.0 
4.0 


-564.0 

-560.8 

556.0 


533.6 
530.6 
526.0 








T/K-- 


= 224.5520 










11 
12 
13 


1473.3901 
3285.0222 
5003.5882 


1473.7176 
3285.4294 
5004.1509 


1473.4993 
3285.1579 
5003.7758 


0.1404 
0.1608 
0.1922 


0.1244 
0.1441 
0.1724 


10.8 
5.1 
4.0 


-563.2 
-560.1 
-555.0 


532.9 
530.0 
525.1 








r/K = 


= 238.6638 










11 
12 
13 


1473.0799 
3284.3317 
5002.5368 


1473.4069 
3284.7381 
5003.0984 


1473.1889 
3284.4672 
5002.7240 


0.1415 
0.1618 
0.1936 


0.1256 
0.1455 
0.1739 


10.8 
5.0 
3.9 


-562.5 
-559.1 
-554.0 


532.2 
529.0 
524.2 








T/K = 


= 252.9801 










11 
12 
13 


1472.7597 
3283.6186 
5001.4510 


1473.0862 
3284.0246 
5002.0119 


1472.8685 
3283.7539 
5001.6380 


0.1425 
0.1632 
0.1952 


0.1267 
0.1468 
0.1755 


10.7 
5.0 
3.9 


-561.7 
-558.7 
-553.5 


531.5 
528.6 
523.7 








T/K-- 


= 273.1600 










11 
12 
13 


1472.2981 
3282.5908 
4999.8863 


1472.6239 
3282.9959 
5000.4459 


1472.4067 
3282.7258 
5000.0728 


0.1440 
0.1648 
0.1971 


0.1283 
0.1486 
0.1777 


10.7 
4.9 
3.9 


-560.7 
-557.6 
-552.3 


530.5 
527.6 
522.6 








T/K = 


= 287.4650 










11 
12 
13 


1471.9652 
3281.8494 
4998.7574 


1472.2905 
3282.2537 
4999.3161 


1472.0736 
3281.9842 
4998.9436 


0.1450 
0.1659 
0.1984 


0.1294 
0.1499 
0.1792 


10.6 
4.9 
3.8 


-560.0 
-556.7 
-551.6 


529.8 
526.7 
521.9 








T/K-- 


= 302.9301 










11 
12 
13 


1471.6001 
3281.0362 
4997.5193 


1471.9249 
3281.4399 
4998.0772 


1471.7084 
3281.1708 
4997.7053 


0.1461 
0.1672 
0.1998 


0.1306 
0.1513 
0.1809 


10.5 
4.9 
3.8 


-559.2 
-556.0 
-550.9 


529.1 
526.1 
521.3 
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Table A2. Mean values and relative standard deviations of (u/ao) for 
argon from 5 modes 



Table A2. Mean values and relative standard deviations of (u/ao) for 
argon from 5 modes — Continued 



p/kPa 


{u/ao)/s-' 


lOV 


p/kPa 


{u/ao)/s-' 


lOV 


7= 302.9166 K 


uo = 88.94344 mm 




7=253.1500 K 


ao = 88.87387 mm 




556.997 


3649.8707 


0.9 


470.069 


3334.3519 


0.4 


556.997 


3649.8707 


0.9 


422.800 


3334.3271 


0.4 


499.559 


3649.3006 


0.8 


375.753 


3334.3074 


0.3 


450.125 


3648.8134 


0.8 


328.542 


3334.2965 


0.3 


400.000 


3648.3241 


0.7 


282.174 


3334.2842 


0.3 


349.777 


3647.8375 


0.6 


234.896 


3334.2783 


0.4 


299.646 


3647.3560 


0.6 


234.761 


3334.2796 


0.4 


299.646 


3647.3562 


0.6 


188.108 


3334.2787 


0.3 


249.957 


3646.8834 


0.6 


140.855 


3334.2838 


0.3 


199.729 


3646.4091 


0.7 


105.364 


3334.2902 


0.4 


149.930 


3645.9430 


0.6 


71.751 


3334.2987 


0.3 


149.930 


3645.9440 


0.5 


49.850 


3334.3070 


0.4 


104.875 


3645.5260 


0.3 


24.921 


3334.3162 


0.4 


74.848 


3645.2492 


0.6 








49.898 


3645.0222 


0.4 


7= 234.3156 K 


ao = 88.84858 mm 




24.974 


3644.7939 


0.5 














480.511 
480.511 


3206.1452 
3206.1468 


0.4 


7-= 293.1300 K 


ao = 88.92949 mm 




0.4 








450.963 
400.725 


3206.2935 
3206.5511 


0.6 


537.608 


3590.1231 


0.6 


0.2 


483.268 


3589.6591 


0.6 


350.927 


3206.8116 


0.7 


430.165 


3589.2111 


0.4 


299.286 


3207.0896 


0.4 


375.878 


3588.7586 


0.5 


250.141 


3207.3522 


0.6 


322.145 


3588.3150 


0.5 


200.101 


3207.6324 


0.4 


267.906 


3587.8722 


0.6 


150.081 


3207.9147 


0.2 


267.906 


3587.8725 


0.7 


100.575 


3208.2014 


0.6 


213.940 


3587.4378 


0.7 


74.711 


3208.3512 


0.6 


159.913 


3587.0067 


0.5 


49.933 


3208.4986 


0.5 


105.775 


3586.5809 


0.3 


25.010 


3208.6469 


0.6 


105.775 


3586.5811 
3586.3436 


0.3 
0.5 








75.145 


7-= 217.0950 K 


ao = 88.82605 mm 




49.763 


3586.1472 
3585.9553 


0.5 
0.4 








24.944 


357.289 


3085.0067 


0.3 








326.549 
291.181 


3085.3776 
3085.8025 


0.5 


r= 273.1600 K 


00 = 88.90141 mm 




0.5 








255.374 
220.638 


3086.2390 
3086.6611 


0.5 


500.233 


3464.8475 


1.1 


0.5 


451.426 


3464.5995 


0.9 


185.789 


3087.0880 


0.6 


399.112 


3464.3411 


0.7 


149.955 


3087.5313 


0.6 


349.798 


3464.1033 


0.7 


149.955 


3087.5317 


0.8 


300.659 


3463.8680 


0.6 


115.584 


3087.9581 


0.5 


250.247 


3463.6343 


0.7 


102.231 


3088.1273 


0.6 


200.970 


3463.4087 


0.4 


102.144 


3088.1284 


0.6 


150.375 


3463.1832 


0.5 


74.809 


3088.4694 


0.3 


127.709 


3463.0844 


0.5 


50.272 


3088.7773 


0.3 


100.288 


3462.9658 


0.4 


24.969 


3089.0984 


0.3 


100.288 


3462.9657 


0.5 








74.917 


3462.8563 


0.3 








50.006 


3462.7512 


0.6 








25.030 


3462.6469 


0.4 
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Table A3. Mean values and standard deviations of (u/uo) for xenon 
corrected for the virtual leak 



p/kPa 


{u/ao)/s-' 


lOV 


7= 302.9166 K 


uo = 88.94344 mm 




302.569 


1998.9335 


0.5 


269.851 


2000.1990 


0.6 


239.836 


2001.3569 


0.7 


199.948 


2002.8931 


0.7 


150.100 


2004.8092 


0.9 


133.515 


2005.4456 


1.0 


105.227 


2006.5302 


0.8 


79.949 


2007.4980 


0.8 


54.958 


2008.4542 


0.7 


44.964 


2008.8364 


0.6 


33.065 


2009.2912 


0.4 


7=273.1600 K 


ao = 88.90141 mm 




299.353 


1893.8741 


0.6 


299.118 


1893.8870 


0.6 


270.053 


1895.4906 


0.7 


239.722 


1897.1588 


0.7 


239.531 


1897.1690 


0.7 


210.107 


1898.7816 


0.7 


179.899 


1900.4307 


0.8 


150.130 


1902.0510 


0.8 


119.991 


1903.6854 


0.8 


89.940 


1905.3097 


0.7 


60.019 


1906.9218 


0.6 


45.059 


1907.7253 


0.5 


29.906 


1908.5384 


0.4 



Acknowledgments 

The full cooperation of the Thermometry Group of 
NIST is gratefully acknowledged. Among many favors 
that we recall, B. W. Mangum checked our gallium 
point cell, G. F. Strouse checked the mercury point cell 
and provided timely calibrations of our thermometers, 
G. T. Furukawa advised us on the preparation of the 
mercury point cell, and D. C. Ripple read this 
manuscript with great care. J. W. Schmidt of the Pres- 
sure and Vacuum Group calibrated our pressure gauge. 
The participation of C. W. M. in this work was possible 
because of the award of an NRC/NBS Postdoctoral 
Research Associateship. The participation of S. J. B. in 
this work was possible because of the award of a NATO 
fellowship by the Science and Engineering Research 
Council (U.K.). This work builds on the efforts of oth- 
ers at NBS/NIST and visitors who contributed to devel- 
opment of the gas constant apparatus itself and to 
understanding of its properties. We gratefully acknowl- 
edge the efforts of M. Waxman, M. J. Greenspan, J. B. 
Mehl, R. S. Davis, T. J. Edwards, J. P. M. Trusler, and 
M. B. Ewing. 



12. References 

[1] M. R. Moldover, M. Waxman, and M. Greenspan, Spherical 
acoustic resonators for temperature and thermophysical prop- 
erty measurements, High Temp. -High Pres. 11, 75-86 (1979). 
[2] J. B. Mehl and M. R. Moldover, Precision acoustic measure- 
ments with a spherical resonator: Ar and C2H4, J. Chem. Phys. 
74 (7), 4062-4077 (1981). 
[3] M. R. Moldover, J. B. Mehl, and M. Greenspan, Gas-filled 
spherical resonator: Theory and experiment, J. Acoust. Soc. 
Am. 79 (2), 253-276 (1986). 
[4] M. R. Moldover, J. P. M. Trusler, T. J. Edwards, T. J. Mehl, and 
R. S. Davis, Measurement of the universal gas constant R using 
a spherical acoustic resonator, J. Res. Natl. Bur. Stand. (U.S.) 93 
(3), 85-114(1988). 
[5] J. B. Mehl and M. R. Moldover, Measurement of the ratio of the 
speed of sound to the speed of light, Phys. Rev. A34 (4), 3341- 
3344 (1986). 
[6] M. B. Ewing, J. B. Mehl, M. R. Moldover, and J. P. M. Trusler, 
Microwave measurements of the thermal expansion of a spher- 
ical cavity, Metrologia, 25 (4), 211-219 (1988). 
[7] M. R. Moldover and J. P. M. Trusler, Accurate acoustic ther- 
mometry I: the triple point of gallium, Metrologia 25, 165-187 
(1988). 
[8] C. W. Meyer and M. R. Moldover, Notes for CODATA meeting 

(impublished). 
[9] A. R. H. Goodwin, Thermophysical Properties from the Speed 
of Soimd, Ph. D. Thesis, University College London (1988) pp. 
113-117. 

[10] M. B. Ewing, M. B. McGlashan, and J. P. M. Trusler, The 
temperature-jump effect and the theory of the thermal boundary 
layer for a spherical resonator. Speeds of sound in argon at 
273.16 K, Metrologia 22, 93-102 (1986). 

[11] M. B. Ewing and J. P. M. Trusler, Primary Acoustic Thermome- 
try Between 100 K and 300 K. Report to the 18th Meeting of the 
CCT, paper CCT/93-28, 20 August 1993. 

[12] R. C. Kemp, W. R. G. Kemp, and L. M. Besley, Determination 
of Thermodynamic Temperatures and Measurements of the Sec- 
ond Virial Coefficient of ''He Between 13.81 K and 287 K using 
a Constant- Volume Gas Thermometer, Metrologia 23, 61-86 
(1986/87). 

[13] D. N. Astrov, L. B. Beliansky, Y. A. Dedikov, S. P. Polunin, and 
A. A. Zakharov, Precision Gas Thermometry between 2.5 K and 
308 K, Metrologia 26, 151-166 (1989); D. N. Astrov, L. B. 
Belyansky, and Y. A. Dedikov, Correction of the gas-ther- 
mometry scale of the VNIIFTRI in the range 2.5 K to 308 K, 
Metrologia 32, 393-395 (1995/96). 

[14] J. E. Martin, T. J. Quinn, and B. Chu, Further measurements of 
thermodynamic temperature using a total radiation thermome- 
ter: the range -130 °C to -1-60 °C, Metrologia 25, 107-112 
(1988). 

[15] L. A. Guildner and R. E. Edsinger, Deviation of the International 
Practical Temperatures from Thermodynamic Temperatures in 
the Temperature Range from 273.16 K to 730 K, J. Res. Natl. 
Bur. Stand. (U.S.) Sec. A 80, 703-737 (1976). 

[16] James B. Mehl, Resonance frequencies of deformed spherical 
resonators, J. Acoust. Soc. Am. 71, 1109-1113 (1982). 

[17] James B. Mehl, Acoustic resonance frequencies of deformed 
spherical resonators 11, J. Acoust. Soc. Am. 79, 278-285 (1986). 

[18] A. F. Clark, G. E. Childs, and G. H. Wallace, Low-temperature 
electrical resistivity of some engineering alloys, in Advances in 
Cryogenic Engineering, Vol. 15, K. D. Timmerhaus, ed.. 
Plenum Press, New York and London (1970) pp. 85-90. 



45 



Volume 104, Number 1, January-February 1999 

Journal of Research of the National Institute of Standards and Technology 



[19] Donald Peckner and I. M. Bernstein, Handbook of Stainless 
Steels, McGraw-Hill Book Company, New York (1977) p. 19-3, 
p. 12-33, and Fig. 46 on p. 12-37. 

[20] H. Preston-Thomas, P. Bloembergen, and T. J. Quinn, Supple- 
mentary information for the international temperature scale of 
1990, Bureau International des Poids et Mesures, December 
1990. 

[21] D. R. White, The contribution of uncertainties in resistance 
measurements to uncertainties in the realisation of ITS-90, 
Preprint for the proceedings of the EUROMET workshop. Un- 
certainties in Temperature Measurements, Paris, 23-24 March 
1998. 

[22] J. L. Riddle, G. T. Furukawa, and H. H. Plumb, Platinum Resis- 
tance Thermometry, NBS Monograph 126, U.S. Govt. Printing 
Office, Washington, DC (1973). 

[23] G. T. Furukawa, J. L. Riddle, W R. Bigge, and E. R. Pfeiffer, 
Application of some metal SRM's as thermometric fixed points, 
Natl. Bur. Standards (U. S.) Special Publication 260-77, U.S. 
Govt. Printing Office, Washington, DC (1982). 

[24] R. A. Aziz, and M. J. Slaman, The argon and krypton inter- 
atomic potentials revisited. Molecular Physics 58, 679-697 
(1986). 

[25] G. C. Maitland, M. Rigby E. B. Smith, and W A. Wakeham, 
Intermolecular Forces, Clarendon Press, Oxford (1981) p. 571. 

[26] A. K. Dham, W J. Meath, A. R. Allnatt, R. A. Aziz, and M. J. 
Slaman, XC and HFD-B potential energy curves for Xe-Xe and 
related physical properties, Chem. Phys. 142, 173-189 (1990). 

[27] R. A. Aziz, private communication. 

[28] M. B. Ewing and J. P. M. Trusler, Interaction second acoustic 
virial coefficients of (N2 + Ar) between 90 and 373 K, Physica 
A184, 437^50 (1992). 

[29] J. H. Dymond and E. B. Smith, The Virial Coefficients of Pure 
Gases and Mixtures: A Critical Compilation, Clarendon Press, 
Oxford (1980) pp. 1-8. 

[30] A. Michels, T. Wassenaar, and P. Louwerse, Isotherms of xenon 
at temperatures between "C and 150 °C and at densities up to 
515 amagats (pressures up to 2800 atmospheres), Physica 20, 99 
(1954). 

[31] H. M. Ledbetter, W F Weston, and E. R. Naimon, Low-temper- 
ature elastic properties of four austenitic stainless steels, J. Appl. 
Phys. 46, 3855-3860 (1975). 

[32] P. R. Bevington, Data Reduction and Error Analysis for the 
Physical Sciences, McGraw-Hill Book Co., New York (1969). 

[33] S. J. Boyes, The speed of soimd in gases with application to 
equations of state and sonic nozzles. Ph. D. Thesis, University 
College London (1992). 

[34] M. B. Ewing and J. P. M. Trusler, Second acoustic virial coeffi- 
cients of nitrogen between 80 and 373 K, Physica A184, 415- 
436 (1992). 

[35] S. J. Boyes, unpublished. 

[36] M. B. Ewing, A. A. Owusu, and J. P. M. Trusler, Second acoustic 
virial coefficients of argon between 100 and 304 K, Physica 
A156, 899-908 (1989). 

[37] M. B. Ewing and A. R. H. Goodwin, An apparatus based on a 
spherical resonator for measuring speed of sound in gases at high 
pressures. Results for argon at temperatures between 255 K and 
300 K and at pressures up to 7 MPa, J. Chem. Thermodynamics 
24, 531-547 (1992). 

[38] A. R. H. Goodwin and M. R. Moldover, Thermophysical proper- 
ties of gaseous refrigerants from speed of sound measurements. 
I. Apparatus, model, and results for 1,1,1,2-tetrafluoroethane 
R134a, J. Chem. Phys. 93, 2741-2753 (1990). 

[39] N. E. Holden, R. L. Martin, and I. L. Barnes, Isotopic composi- 
tions of the elements 1983, Pure Appl. Chem. 56, 675-694 
(1984). 



About the authors: Dr. Michael R. Moldover is a NIST 
Fellow and is also the Leader of the Fluid Science 
Group in the Physical and Chemical Properties Divi- 
sion of the NIST Chemical Science and Technology 
Laboratory. He has worked on a wide range of thermo- 
physical properties measurements including a redeter- 
mination of the universal gas constant R using acoustic 
methods. 

Dr Steven J. Boyes was an SERC/NATO Postdoctoral 
Fellow when this research was conducted. He is now a 
physical chemist working on noncontact thermometry 
in the Centre for Basic and Thermal Metrology of the 
National Physical Laboratory. 

Dr. Christopher W. Meyer was an NRC Postdoctoral 
Fellow when this research was conducted. He is now a 
physicist in the Thermometry Group of the Process 
Measurements Division in the NIST Chemical Science 
and Technology Laboratory and works on realizing 
ITS-90 at low temperatures and on temperature mea- 
surements during rapid thermal processing. 

Dr. Anthony R. H. Goodwin was a physical chemist 
in the Fluid Science Group when this research was 
conducted. He is now a Research Scientist at Schlum- 
berger Cambridge Research. 



46 



